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f) to the 2007 world np data below 350 MeV. The wave functions are expanded in a series of 
analytical functions (with the correct asymptotic behavior at both large and small arguments) that 
can be Fourier-transformed from momentum to coordinate space and are convenient to use in any 
application. A fortran subroutine to compute these wave functions can be obtained from the authors. 
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I. INTRODUCTION 

This paper presents the effective range expansions and 
relativistic deuteron wave functions that accompany the 
recent high precision fits to the 2007 world np data (con- 
taining 3788 data) below 350 MeV [11. These fits, ob- 
tained using the covariant spectator theory (GST) OH], 
resulted in two new one boson exchange models of the np 
interaction, both of which also reproduce the experimen- 
tal triton binding energy without introducing additional 
irreducible three-nucleon forces ^. One model (WJC-1) 
has 27 parameters and fits with a x^/iVdata = 1.06. The 
other model (WJC-2) has only 15 parameters and fits 



with a X /N, 



data 



1.12. 



The main body of the paper is divided into three sec- 
tions. In Sec. II the 6th order effective range expansions 
for the ^Sq and '^S'l phase shifts are given; it is shown 
that these expansions provide an excellent description of 
the phase shifts up to 50 MeV lab kinetic energy. Then, 
in Sec. Ill, the relativistic deuteron wave functions that 
automatically emerge from these fits are discussed. The 
deuteron binding energy was constrained during the fits, 
so both models have the correct binding energy. 

Section III begins with a review of the definitions of 
the wave functions and a description of how they are re- 
lated to the relativistic dnp vertex functions defined in 
the GST. It is shown how the CST wave functions are 
decomposed into four independent amplitudes, and how 
these can be identified with the familiar u (S-state) and 
w (D-state) components, plus two P-state components of 
purely relativistic origin [5l|6]. Because these P-states are 
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associated with the virtual antinucleon degrees of free- 
dom they have positive parity. One P-wave component 
has a symmetric total spin triplet structure (denoted Vt) 
and the other an antisymmetric total spin singlet struc- 
ture (denoted Vs). The antisymmetric Vs wave function 
would be zero if both nucleons had the same energy, but 
because one is on-shell and the other off-shell, the rel- 
ative energy is Ep — \Md (where Ep is the total rel- 
ativistic energy of a nucleon with three momentum p 
and Md is the deuteron mass), allowing for the existence 
of this odd state. The normalization of the wave func- 
tions is discussed [7]. It turns out that the probability 
of the WJC-1 P-state components is only about 0.3% 
while the WJC-2 P-state components are much smaller 
(with a combined probability of less than 0.02%). These 
components are retained because they are required by the 
manifest covariance of the CST, but they are not present 
in some other relativistic approaches [8] . The wave func- 
tions are expanded in a series of analytic functions that 
can be conveniently Fourier-transformed to coordinate 
space, giving analytical wave functions in both momen- 
tum and coordinate space convenient to use in any ap- 
plication. The section concludes with a discussion of the 
asymptotic D/S ratio. Calculations of the magnetic and 
quadrupole moments require the evaluation of interac- 
tion currents, which is beyond the scope of this paper 
and is deferred to future work. 

Our results are summarized in the conclusion. Sec. IV. 
There are also three Appendices. Appendix \K\ gives de- 
tails of the construction and evaluation of the kernels 
omitted from Ref. I [T|. Appendix[B| discusses the method 
used to evaluate the angular integrals, complicated by the 
presence of a cusp in the region of integration. Finally, 
Appendix [C] gives more details about the definitions, ex- 
traction, and normalization of deuteron wave functions. 



II. EFFECTIVE RANGE EXPANSIONS 

The four-term effective range expansion for the ^Sq 
and '^S'l phases that we use is 



k cot (5 = — 



I + ro^ (i + {rokf [P + Q [rokf] ) (2.1) 



where a is the scattering length, tq the effective range, 
and P and Q the dimensionless third and fourth order 
parameters. Here 



k^ = ^mEi 



ab 



(2.2) 



where k is the magnitude of the cm momentum of the nu- 
cleons, -Eiab the laboratory kinetic energy, and we used 
m = 938.9 MeV and hc= 197.3288 MeV fm = 1. RecaU 



that Eq. ( 2.2 1 holds for both relativistic and nonrelativis- 



tic kinematics, allowing us to compare our results directly 
with nonrelativistic calculations. 

To fix these parameters we evaluated the phase shifts 
at four energies: 0.0001, 10, 25, and 50 MeV. The first 
energy of 0.0001 MeV is so low that it effectively fixes 
the scattering length (and hence the total cross section at 
the np threshold) . The use of these four energies resulted 
in a small error between the phases calculated at 1 and 
5 MeV, and those predicted using (2.1). The effective 
range parameters and the errors at 1 and 5 MeV are 
summarized in Table [J 

The effective range expansion is surprisingly accurate 
over a wide energy range. As shown in Fig. [T] the expan- 
sions are qualitatively accurate up to E^iab of about 150 
MeV (the small errors shown in Table IT] are completely 
invisible on the figure). In the fitting we used the expan- 
sions instead of the phase shifts to calculate the S wave 
contributions to the observables of all energies below 50 
MeV. In doing the final minimizations we constrained 
both the deuteron binding energy and the ^ Sq phase shift 
at 0.0001 MeV (which we fixed at S = 1.4937°). This was 
a very effective way to maintain an accurate fit to the low 
energy cross section data. 
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FIG. 1. (Color online) The effective range expansions for the 
S-waves (smooth line) compared to calculated phase shifts (solid 
dots). The phases are in degrees. 



TABLE I. The effective range parameters defined in 
Eq. ( |2.1[ ) and the fractional errors A at 1 and 5 MeV 
for each model. The fractional error is defined to be 
A = (cot (5 — cot (5cfl) /cot (5, where ScS is the phase 
predicted by Eq. (2.1 1 and S the phase emerging from 



solutions of the scattering equation. 



'So 



^5i 



parameter 



WJC-1 WJC-2 WJC-1 WJC-2 



a (fn 



-23.7494 -23.7496 5.4295 5.4342 



ro (fm) 


2.6261 


2.6623 


1.7601 


1.7666 


P 


0.0007 


-0.0005 


0.0172 


0.0206 


Q 


0.0033 


0.0032 


0.0060 


0.0058 


A (1 MeV) 


0.0044 


0.0038 


-0.0011 


-0.0009 


A (5 MeV) 


0.0033 


0.0030 


-0.0012 


-0.0012 



III. DEUTERON WAVE FUNCTIONS 



The CST deuteron wave functions have been devel- 
oped and defined in a number of references. The earliest 
reference is the work of Buck and Gross [5]; more re- 
cent references include Refs. [9] (referred to as Ref. II) 
and Ref. [10], where many practical details for how to 
use the wave functions in practical calculations are de- 
veloped. While all of these references use the same basic 
definitions, details of how to interpret and use the wave 
functions have improved with time. In section we first 
present a brief review of the necessary definitions (relat- 
ing them to the earlier references), and then present the 
numerical results for the wave functions of models WJC-1 
and WJC-2. 



A. Definitions of the vertex function and tiie 
covariant wave function 

The equation satisfied by the np scattering amphtude 
is (in the notation of Eq. (2.3) of Ref. I) 

Mi2(p,p';F)=Fi2(p,p';P) 

Fi2(p, k; P)G2{k2)Mi2{k,p'; P) , (3.1) 

where P is the conserved total four- momentum, and p, p', 
and k are relative four-momenta related to the momenta 
of particles 1 and 2 by pi = |P + p, p2 = \P — P- 
Ek^ = \J°m?- -I- k^ is the energy of the on-shell particle 1, 
and the covariant integral is 



d^k 



(27r)3 Ek 



(3.2) 



Note that these covariant operators can be written ei- 
ther in terms of the independent momenta {P,p,p'} or 
{P,pi,p']^}. The scattering amplitude, 

Mi2(p,p';P) = iV4^v„/3/^'(P,p';-P) 

= MQ(pi,Ai)A^aQ',0/3'(p,p';P)"Q'(pi,Ai), (3.3) 

is the matrix element of the Feynman scattering ampli- 
tude M. between positive-energy Dirac spinors of particle 
1. The definition of the nucleon spinors u(pi, Ai) (with 
Ai the helicity of the nucleon) is given in Eq. ( |A3[ ); the 
on-shell spinor has four-momentum px with three mo- 
mentum component pi, and p\ — \ 
for the ofF-shell particle 2 is 



The propagator 



(m + 1^2) 



^2 



ie 



h\k2) (3.4) 



G2(fc2)^G,,,(fc2) 

with fc2 — P — ki, kf = m? , and h{k2) [denoted H{k2) in 
Ref. I] the form factor of the off-shell nucleon (related to 
its self energy), normalized to unity when fc| = m? 



h{k2) = 



{A%-my 



-\ 2 



(A2^ - m2)2 + (m2 - A:2)2 



(3.5) 



See AppendixlOfor further discussion of the nucleon form 
factor h and the role it plays in the deuteron wave func- 
tions. 

Just below threshold, the np scattering amplitude has 
a pole at the deuteron mass, P^ = Mj. Near the pole 
this amplitude, in the CST, can be written 



Mx,x'^.jip,ip,p';P) = Rx,x[,p0'{p,p';P) 

^gx,fi{p,P,\d)Gx[p'{p\PAd) 



Ad 



M! 



p2 



(3.6) 



where Ai, X'l are the helicities of the outgoing and incom- 
ing on-shell particle 1, /3,/3' are the Dirac indices of the 
off-shell particle 2, A^ is the helicity of the deuteron, Q is 




FIG. 2. (Color online) Diagrammatic representation of the vertex 
function ||3.8|l. The on-shell particle is labeled by the X, so that 



Pj = m^, and the relative four-momentum is p = , (Pi ~ Pa). The 



deuteron is also on-shell, with P^ 



■Ml 



the contracted dnp vertex functions describing the cou- 
pling of the deuteron to the neutron and proton, and R is 
a remainder function finite at the pole. The equation for 
g can be found by substituting ( |3.6[ ) into the scattering 
equation ( |3.1[ ) and demanding that it hold at the pole 
(see Refs. [3] [7] and II), giving 

^Ai,9(P,P,Ad) =- / Vx,X[,l3aiP,k]P) 
Jk[ 

^G^p,{k2)gx[p'{k,P,Xd). (3.7) 

The contracted vertex function is a Lorentz scalar 
product of a Dirac matrix element and the covariant po- 
larization vector of the deuteron, ^^(A;;): 



gxAp,p,Xd)^oi^ip,p)U>^^), 



(3.8) 



where O is the "uncontracted" dnp vertex function for 
particle 2 off-sheU [c.f. Eq. (2.27) of Ref. II] 



OlJp,P)^ T^ip,P)C 



pp 



4(pi,Ai), (3.9) 



with C the Dirac charge conjugation matrix. To simplify 
the language, g will always be called the contracted ver- 
tex function, and the functions O (or sometimes F) will 
be referred to simply as the "vertex function" . The con- 
tracted vertex function is illustrated diagramatically in 
Fig.|2l 

The relativistic deuteron wave functions were defined 
previously in Eq. (2.32) of Ref. II and Eq. (2.29) of 
Ref. [10]: 

i^x,PM{P,P) = NdGiip,{p2)gx,fi'{p,P,Xd). (3.10) 

Here N^ is a normalization constant, chosen to be 

1 



N, 



V(2^)3 2Md 



(3.11) 



As discussed in Ref. [TD], it is the vertex function that 
enters directly into any Feynman diagram involving an 
incoming or outgoing deuteron state, so if a Feynman 
amplitude is expressed in terms of the wave function in- 
stead of the vertex function, one must be careful to divide 



by the normalization constant (3.11 1. The reason for this 



choice of normalization constant will be discussed further 
below. 



The normalization of the wave function (3.10) is dis- 
cussed in detail in Refs. OIT], II, and Appendix [C| The 
exact result, in the notation of Eq. (3.7), is 



1 



/,t 



XV'A,0,A'(P',P)]>, (3.12) 



d ^ -'Pi 
Pi "'P'l 



where summation over repeated indices is implied, and 

AFAiAi,/j/3'(p,p';^) = ^Vx,x[,pf,'{p,p;P), (3.13) 

with the partial derivative holding pi and p[ constant. 
This normalization integral will be simplified and re- 
lated to the usual nonrelativistic normalization in Ap- 
pendix |C3| below. 

While it is often convenient to use the fully covariant 



expression (3.10), the physical content of the wave func- 



tion can be displayed if we expand the off-shell propaga- 
tor into positive and negative energy pieces. (Note that 
an alternative expansion given in Eq. (A21) of Ref. [TU] 
is sometimes more convenient for applications.) Working 
in the rest frame, where we choose pi — — P2 — p, the 
propagator can be decomposed 



GapiP2) 



E, 



E 



Ua("P, A2)u/3(-p, A2) 



i;a(p, A2)t;/3(p, A2) 



Er 



Y, G"^ ip)< (-P, A2)S^= (-P, A2) (3.14) 



A2P2 



where p2 is the p-spin of the off-shell particle 2. The 
p-spin propagators G''{p) = G'^{p) are 



G+(p) 



1 



2E„ - M, 



jp 



d 



G-{P) = - 



1 



(3.15) 



and u and v are the positive- and negative-energy spinors 
of particle 2, with u^ — u^ and m+(— p, A2 ) = m(— p, A2) 
and u~(— p, A2) = w(p,A2) [see Eq. (A3)]. Using this 



notation the wave function can be conveniently written 

V'Aia,A,(p,P) - Y. Ca„A.(pX(-P'A2), (3.16) 

where the wave functions -0''^ are 

^fiA.,A.(p) = N,^G''^{p)Ti\^^^,^{^). (3.17) 

Suppressing the Dirac indices in the matrix elements 
gives 

rAiA..A.(P)-"^^(-P>A2)t;Ai/3(p,P,A,) 

= S^H-p, A2)r^(p, P)CS^(p, Ai) U^d). (3.18) 



These definitions agree with Eq. (2.67) of Ref. II. Note 
that r^^_^ y^ (p) contains no additional factors of m/E; 
these are written explicitly in (3.17), and will be identi- 
fied whenever they appear below. 



B. 



Partial wave CST equations for the deuteron 
vertex functions 



Equation (3.1), and the companion equation (3.7) for 



the bound state are solved by expanding the amplitudes 
into partial waves. The connection between the wave 
functions in momentum space and their partial waves will 
be briefly outlined in this subsection; for further details 
see Appendix E of Ref. I, and, for the deuteron channel. 
Appendix [^ 

The first step is to extract the J = 1 partial waves from 
the momentum space amplitudes F^^^ ^ (p) defined in 



Eq. (3.18). These amplitudes have a simple angular de- 
pendence. If p = p{sin0cos0, sin6'sin(/), cos6'} (where, 
from here up to section |III G| below and again in Ap- 
pendix O we will use the notation p = |p| in order to 
avoid confusion with the four-momentum p) , the form of 
these functions is 



'^x',Xn,\AP) 



Di:,{^,e,0)Til,^^^^{p), (3.19) 



where A = Ai — A2, A^ = Mj is the polarization of the 
deuteron, and D^ are the J = 1 rotation matrices. (This 
can be obtained from Eq. (E29) of Ref. I by recalling that 
the scattering amplitude at the deuteron pole is propor- 
tional to the product of vertex functions for the incoming 
and outgoing nucleons.) In order to keep our notation 
succinct, we suppress the usual index J on the partial- 
wave vertex function T^ ''^ ^ (p) . An argument p instead 
of p identifies it unambiguously as a partial-wave vertex 
function . If w e choose p to lie in the xz plane (cj) — 0), we 
19 ) to find the connection between T'^ ^ ^ (p) 
ciiiu 1- X X \ (p) . Using the normalization condition for 
the d""s 



1 



i(2J+l)/ sm9de[di,x{0)] , (3.20) 



the relation is 



pP2 

Ai A2 , A 



(p) = VSn / sin6' 
' Jo 



d0 



xdL,A(e)rf,,,,,,(p) 



(3.21) 



and the deuteron wave functions corresponding to (3.17) 
are 

TO 



Ca2.a.(p) 



N, 



d — GP^p)TP\ . 



(p). (3.22) 



The next step is to use the properties of the partial 
wave amplitudes under parity and particle interchange 
(derived in Refs. I and II) to show that there are only 
four independent vertex functions. We begin by restor- 
ing reference to the p-spin of particle 1 and the relative 



energy po = (pio^P2o)/2 (previously suppressed for con- 
venience) , which transforms the partial wave vertex func- 
tions (3.21 1 into 



pP2 



,ip)^^tX,xMPo)^ (3-23) 

where, in the cm frame with particle 1 on shell, po — 
Ep — Md/2 > 0. Under parity (V) and particle inter- 
change (7^12 ), it was shown in Appendix E of Ref. I that 
these partial wave amplitudes (for J = 1) satisfy the 
transformation properties 

A,(P.Po) = P2r+^'_A2,A,(P'Po) 

(p,Po)=rf+.Ajp,-Po). (3.24) 



-pnx 



r+P2 
2-^ AiA2,Ad 



The parity relation shows that only two (of the four pos- 
sible) helicity states are independent (which we chose to 
be Ai = +; the states with Ai = — can be obtained 
using parity) and the exchange relation shows that the 
exchange amplitudes (with particle 2 on-shell and neg- 
ative relative energy) can be obtained from the direct 
amplitudes (with particle 1 on-shell and positive rela- 
tive energy). The detailed construction of the deuteron 
bound state equation is reviewed in Appendix |C 1[ where 
it is shown how these symmetry properties are derived 
and used to reduce the initial set of coupled equations to 
only four independent ones. 



Using these results, the vertex functions (3.21 ) can be 
organized into the column vector 

r;+,A,(p) 

r:-,A.(p) 
Vr;_,.,(p)y' 

where, in the matrix representations, we return to the 
notation p to represent the momentum dependence of the 



\Tx,ip))=N, 



E„ 



(3.25) 



vertex function. With this notation, the original Eq. (3.7) 



can now be written as four independent equations for 
these vertex functions. In a convenient matrix form the 
equations are 



\TxAp)) 



edkYd{p,k)g{k)\TxAk)), (3.26) 



where Vd(p, k) (the kernel at the deuteron pole) and g{k) 
are 4x4 matrices. The equations are independent of the 
projection of the total angular momentum, which for the 
deuteron is its helicity {mj = Ad). With the factors of 
m/E explicitly included in the definition of the column 
vector (3.25) and in the kernel Vd(p, k), the propagator 



is the 4x4 diagonal matrix: 
/ G+{k) 



g(fc) 







\ 



G-{k) 

G+{k) 

V G'{k) J 



(3.27) 



with Gf{k) defined in Eq. (3.15), and the kernel is the 
4x4 matrix 



/«9 



Vd(p,fc) = 



\ 



,++ 



\v 



16 



(3.28) 



/ 



where the wf^''^ were previously defined in Eq. (E45) of 
Ref. I and also in Eq. (C6). 



C. Decomposition into the S, D, and P states 



The Dirac operator F'' [from Eq. (3.9)] can be ex- 
panded in terms of four scalar invariant functions, 
F, G, H, and / as given in Eq. (33) of Ref. ^ and 
Eq. (2.28) of Ref. II, 



r^(p,p) ^FY 



--pf^ 

m 
(to— ^) 



{hY 



-p' 



(3.29) 



These functions differ by some factors of 2 from those 
originally introduced by Blankenbecler and Cook [6]. 
Since the deuteron is on mass-shell, so that P^ = W"^ = 
Af J , they are functions of p^ only. 

Using the explicit definition of the spinors given in Ap- 
pendix ffl and substituting the expansion (3.29) into the 
Dirac matrix elements (3.18), the wave functions (3.17) 



can be reduced to 2-component matrix elements, with 
four independent scalar wave functions u,w,Vt, and Vs 
written as linear combinations of the four invariant func- 
tions F,G,iJ, and/ [5]: 



^AiA.,A.(P) = 
, w(p) 

V2 
.(P) = 

72 



1 



w(p)cr-^ 



i, 



An 
3~ 



IC72 

V2 



Xa]^ 



Ai A2,A 



vAp) p-I, 



(t-P^-^,^-P-^,^ 



ia2 

71 



X.,. (3.30) 



For explicit formulae connecting the F, G, H, and / to 
u,w,vt, and Vs, see Eqs. (45) and (46) of Ref. [5] [note 
that the normalization factor (3.11) sets the scale of the 



connection]. Here p = {sin 6*, 0, cos 6*} lies in the x,z 
plane. For a deuteron at rest with its polarization vector 
defined along the z direction, the three-vector compo- 
nents of the polarization are 




e±=T 




(3.31) 



The nucleoli spinors can be either hehcity spinors, or 
spinors with spin projections along the (fixed) z direc- 
tion. In the latter case, we must replace Ai — >■ s\ and 
— A2 — >■ S2 (since, when 6* = 0, the helicity and spin pro- 
jection of particle 2 are in opposite directions). In either 
case (by inspection) the wave functions are real. 

The wave functions u,w,vt., and Vs have a convenient 
physical interpretation that follows from the structure 
of the 2-component matrix elements that multiply them. 
This follows most directly if the nucleon spin projections 
arc fixed in the z direction. In this case the p — + wave 
functions is symmetric in the spins, 



(<..,A.(P))^ 



<..A.(P) = <.,,A.(P), (3-32) 



showing that the spin state is a spin-one triplet. The 
p = — term includes a symmetric piece (the spin triplet 
vt) and an antisymmetric piece (corresponding to a spin 
singlet state, denoted by the s subscript on Vs). 

In the CST, wave functions in coordinate space are 
defined to be the Fourier transforms of the momentum 
space wave functions (3.301 [5_ 



^siS2,A<iW (27r)3/2 



rf'pe^^"<..,A.(p)- (3-33) 



Using the expansion of the plane wave 



e'P- = 47r^ J2 i'MpT)Y;^{p)YUi) (3-34) 



e m=- 



and noting that the coefficient of w(p) involves a linear 
combination of the Y2m spherical harmonics 



3p^p^-S^^=Y,C%Y2,n{p) 



(3.35) 



(where the Cij are constants) it follows that the coor- 
dinate space spin representation of the deuteron wave 
functions is 



"^^.^.A.W 



vtir) 




(3.36) 



where, denoting the typical wave function by zg (so that 
zq = u, Z2 = w, and zi = vt or w^), the momentum and 
position space wave functions are related by the spherical 
Bessel transforms 



zi{p) 
ztX^) 



rdrj£(pr)z^(r) 



P^dpj({pr)zi{p) (3.37) 



where ji is the spherical Bessel function of order £ with 
the convenient recursion relation 



h{z) 



1 d \ sin z 
z dz ) z 



(3.38) 



Note the appearance of the minus sign multiplying w{v) 
and the factor of i multiplying both of the P-state terms. 
These come from the factor of i^ in the plane wave ex- 
pansion and can be easily overlooked. 



D. Normalization of the wave functions 



In Appendix C3 the normalization condition (3.12) is 
reduced to the following simple form 



p^dp 



u^{p)+w\p) + v^{p)+v^,{p) 



dV 
dm 



(3.39) 



The derivative term is negative and of the order of a few 
percent (see Table III]). 

The origin of the derivative term can be understood in 
two rather different ways. First, it can be derived from 
the requirement that a nonlinear version of the equation, 
used to derive the unitarity relation for positive ener- 
gies, also hold at the deuteron pole [3l [7]. From this 
point of view, the normalization condition ensures that 
the strength of the pole is not altered by repeated inter- 
actions near the pole and shows that the normalization is 
a consequence of the equation itself. But the normaliza- 
tion condition also follows from the requirement of cur- 
rent conservation, where it is seen to give precisely the 
correct factor to insure that the conservation of charge 
is an automatic consequence of the correct normalization 
of the deuteron wave functions f^ . 

The second interpretation shows, indirectly, that in- 
teraction current contributions to the deuteron charge 
(which give rise to the derivative term) are necessarily 
of the order of a few percent. Hence, the calculation of 
any electromagnetic property of the deuteron (including 
the magnetic and quadrupole moments) that does not 
include interaction currents can be expected ab initio to 
be in error by a few percent. Since the famous discrep- 
ancy in the quadrupole moment is about 5%, it is quite 
possible that a careful CST calculation (including inter- 
action currents) could explain it, but this calculation is, 
unfortunately, beyond the scope of this paper. For this 
reason, we will not report deuteron moments here. 

If these wave functions in are used in any calculation 
that neglects interaction currents, it is probably a better 
approximation to use what we will refer to as the scaled 
normalization condition 



1 



p dp 



u''{p) + w\p)+vl{p)+vl{p) . (3.40) 



This is the normalization condition used in the earliest 
treatments of the CST deuteron wave functions fS] , and 



TABLE II. The deuteron probabilities (in percent) . Both 
the exact and scaled probabilities are shown. 



WJC-1 



WJC-2 



probability 



exact 



scaled 



exact 



scaled 



Ps 

Pd 



97.3876 


92.3330 


95.7607 


93.5985 


7.7452 


7.3432 


6.5301 


6.3827 


0.1180 


0.1119 


0.0103 


0.0101 


0.2234 


0.2118 


0.0090 


0.0088 



Pv, 


0.2234 


0.2118 


0.0090 


0.0088 


j:p 

(V) 


105.4743 
-5.4743 


100.0000 


102.3101 
-2.3101 


100.0000 


total 


100.0000 


— 


100.0000 


— 



is also appropriate for electromagnetic calculations using 
the relativistic impulse approximation (RIA) [l0|,lllj. 



E. Asymptotic behavior of the virave functions 

Using the CST equations for the bound state it is possi- 
ble to predict the asymptotic behavior of the wave func- 
tions. The derivation and discussion are given in Ap- 
pendix (Cl The result is: 



u(p) -> 



1 



p" 



u)(p) - 

Ws(p) 



1 



p" 



(3.41) 



Before the convergence of the angular integrals was im- 
proved (as discussed in Appendix mv the numerical so- 
lutions had power law behaviors which differed from the 
integers predicted in ( 3.41[ ) by about 0.1 to 0.2. With 
the newly converged angular integrals, the actual solu- 
tions exhibit the correct behavior. 

The odd powers for u and w (and even powers for the 
v^s) require special consideration when analytic repre- 
sentations of the wave functions are constructed in the 
following sections. 



F. Numerical results for the wave functions 

The solutions for the wave functions are tabulated 
in Tables |III| and |IV[ with the high momentum behav- 
ior shown graphically in Fig. [3] The figure shows that 
the asymptotic estimates (3.41 1 hold for momenta larger 



than about 10 GeV, and that in this large momentum 
region the small P-state wave functions are larger than 
the dominant S and D-state wave functions. The wave 
functions of model WJC-1 satisfy the approximate rela- 
tions u c^ w, and the D-state wave functions of the two 
models are also comparable, but the S-state wave func- 
tion of model WJC-2 has a zero at large momentum and 
is considerably smaller than its counterpart. The P-state 
wave functions for both models are comparable to within 
a factor of 2, with vt ^ —Vg. 



The probabilities for each of the components of the 
wave function arc summarized in Table ITll This Table re- 



ports both the exact probabilities (3.391 and the renor- 
malized probabilities (3.40). 



G. Fitted wave functions 

It is convenient to fit the deuteron wave functions to 
a series of simple functions that can be analytically in- 
tegrated to obtain the wave functions in r space, and to 
interpolate for any value of the momentum. This subsec- 
tion will describe in detail how this is done. 



1. Scaling the wave functions 

The first step in the fitting process is to scale out the 
rapid dependence of each wave function on the momen- 
tum, so that it may be studied on a linear plot. Con- 
structing these scaling functions took some care, but in 
the end we found that the functions 



Usca.lc{p) 
Wscalcip) 



N 



Vscalcip) = Np 



{mil -f p^)'^/^ 

{mil +P^f'^ 
1 



1 



R 



{al+p^) {mli+p^) 

1 R 

{al+p^) {mli+p^) 
R 



(m^p + p2)5/2 (-„^2^ j^ p2)5/2 



(3.42) 



(where the functional form of Ugcaic is used to scale both 
Vt and Ws, and we return to denoting the magnitude of 
the three-momentum by p [instead of p, as we did in the 
previous sections]) work very well. Note that these scal- 
ing functions have the desired (and observed) p^ behavior 
at small p and the correct asymptotic p dependence for 
all the states. The parameter ao is determined from the 
asymptotic relation 



2\ m? + Oq — 2m = e 



(3.43) 



where e — 2.2246 MeV is the deuteron binding energy. 
This gives 



ao 



ie2 



45.7159 MeV 



(3.44) 



which differs slightly from the nonrelativistic ao — \fmt^ 
as has been emphasized in the literature [8j . This param- 
eter is fixed by the deuteron binding energy, but since the 
P-states vanish outside of the range of the potential, their 
leading mass, ttisq is treated as a free parameter. This 
and the other parameters used in the scale functions are 
given in Tables IV] and |VI[ 




FIG. 3. (Color online) Log-log plot of the momentum space wave functions z^ (p), defined in Eq. | |3.30[ l. (Left panel): WJC-1 wave 
functions tabulated in Table nn] compared to the asymptotic limits given in Eq. ( |3.41[ |. Solid circles: uq; solid diamonds: U12; open 
squares: vt; open circles: Vg. The solid lines are proportional to p~^ and p~'*. (Right panel): WJC-2 wave functions tabulated in Table 
|lV[ The solid lines are proportional to p~^ and p~^, and are the same as those shown in the left panel. 



2. Fitting the momentum space wave functions 

In order to obtain simple and accurate fits to the wave 
functions, we found it convenient to use two types of 
expansion functions. For all but the high momentum 
tail, we use functions that go asymptotically like an even 
power of p for the S and D-states (in practice, p~^), and 
an odd power of p for the P-states ( p~^). These func- 
tions go to zero precisely one power of p faster than the 
observed asymptotic behavior and are very conveniently 
transformed to coordinate space giving a superposition of 
simple exponentials. However, the asymptotic tail can- 
not be accurately described by such functions, and it is 
best to use one special function with the correct asymp- 
totic behavior in order to describe the wave functions at 
very large p. This function has a Fourier transform that 
can be written in terms of the modified Bessel functions 
of the second kind, Kn{z). 

For the first set of functions we choose 



Glip) 



pm\ 



M, 



2ne-e 



TT (rnf +p'^){Mf +p'^)^' ' 



(3.45) 



where the factor of \/2/n is introduced for convenience, 
£ is the angular momentum of the state, n(—2 for the S 
and P states and ng — 3 for the D-state. Note that, near 



p = 0, the expansion functions have the normalization 



GKO) 



M,, 



(3.46) 



The tail wave function, chosen to have the correct fall- 
off as p -^ cx) and the same normalization at p = 0, is 



G^P) 



P 



^^2nH 



(M2-t-p2)«, + i 



(3.47) 



Denoting the typical wave function by zi (as we did 
above) the full momentum space wave functions are ex- 
panded in terms of Gi 



zi{p) 



Y.^lG\{p). 



(3.48) 



As this notation implies, the last function in the sum 
(when i — n) IS the "tail" function (3.47), while the first 
n — 1 are of the type (3.45). Since all of the functions 



have the same normalization at p = 0, the relative size 
of the expansion coefficients is a measure of the relative 
size of each function (at least at small momenta). 

The masses m^ and A/^ (for i = 1 to n— 1) used in each 
of the wave functions ( 3.45 ) were defined by the relations 



nii = ai + (i - l)?7ii 
Mi = rrii + ma: 



(3.49) 



TABLE III. Model WJC-1 momentuin space wave functions z^(p), defined in Eq. ( 3.30 1. The momenta are in 
GeV and the wave functions are in units of (GeV)""^". These wave functions are normahzed to 105.4743%. 
The original value of w at the lowest momentum point (-0.690345E-02), which resulted from the imprecise 
cancellation of two very large terms, can be replaced by the bold face number obtained using p^ scaling. 



0.527313E-03 


0.146227E+03 


0.495991E-03 


-0.623843E-03 


-0.358658E-03 


0.277683E-02 


0.145699E+03 


0.137542E-01 


-0.326777E-02 


-0.192465E-02 


0.681763E-02 


0.142975E+03 


0.813298E-01 


-0.802292E-02 


-0.471289E-02 


0.126402E-01 


0.135554E+03 


0.264888E+00 


-0.148513E-01 


-0.871988E-02 


0.202305E-01 


0.121587E+03 


0.607884E+00 


-0.236810E-01 


-0.139049E-01 


0.295712E-01 


0.101834E+03 


0.108549E+01 


-0.343738E-01 


-0.201942E-01 


0.406428E-01 


0.798054E+02 


0.160146E+01 


-0.466979E-01 


-0.274694E-01 


0.534249E-01 


0.594035E+02 


0.204959E+01 


-0.603055E-01 


-0.355616E-01 


0.678972E-01 


0.428005E+02 


0.236963E+01 


-0.747258E-01 


-0.442522E-01 


0.840413E-01 


0.303155E+02 


0.255133E+01 


-0.893826E-01 


-0.532839E-01 


0.101842E+00 


0.213130E+02 


0.261197E+01 


-0.103638E+00 


-0.623822E-01 


0.121288E+00 


0.149401E+02 


0.257733E+01 


-0.116855E+00 


-0.712774E-01 


0.142375E+00 


0.104515E+02 


0.247282E+01 


-0.128471E+00 


-0.796948E-01 


0.165107E+00 


0.728541E+01 


0.232096E+01 


-0.137990E+00 


-0.876522E-01 


0.189495E+00 


0.504321E+01 


0.213982E+01 


-0.145164E+00 


-0.948577E-01 


0.215565E+00 


0.344843E+01 


0.194353E+01 


-0.149827E+00 


-0.101215E+00 


0.243352E+00 


0.231128E+01 


0.174361E+01 


-0.151964E+00 


-0.106954E+00 


0.272908E+00 


0.149976E+01 


0.154713E+01 


-0.151653E+00 


-0.111870E+00 


0.304301E+00 


0.922698E+00 


0.135976E+01 


-0.149087E+00 


-0.115989E+00 


0.337619E+00 


0.515394E+00 


0.118496E+01 


-0.144410E+00 


-0.119420E+00 


0.372971E+00 


0.231152E+0Q 


0.102511E+01 


-0.137607E+00 


-0.122164E+00 


0.410488E+00 


0.350502E-01 


0.881942E+00 


-0.128414E+00 


-0.123858E+00 


0.450334E+00 


-0.994438E-01 


0.755982E+00 


-0.116562E+00 


-0.124047E+00 


0.492699E+00 


-0.189481E+00 


0.645099E+00 


-0.102598E+00 


-0.121827E+00 


0.537814E+00 


-0.242417E+00 


0.544069E+00 


-0.878148E-01 


-0.116934E+00 


0.585950E+00 


-0.262050E+00 


0.448848E+00 


-0.732091E-01 


-0.109318E+00 


0.637426E+00 


-0.255246E+00 


0.360144E+00 


-0.592921E-01 


-0.995041E-01 


0.692620E+00 


-0.230996E+00 


0.280430E+00 


-0.465163E-01 


-0.878985E-01 


0.751979E+00 


-0.196819E+00 


0.211737E+00 


-0.352080E-01 


-0.752255E-01 


0.816027E+00 


-0.159320E+00 


0.154625E+00 


-0.255953E-01 


-0.621228E-01 


0.885387E+00 


-0.122585E+00 


0.109110E+00 


-0.177552E-01 


-0.494042E-01 


0.960798E+00 


-0.900194E-01 


0.742130E-01 


-0.116628E-01 


-0.377210E-01 


0.104314E+01 


-0.631224E-01 


0.485750E-01 


-0.718102E-02 


-0.276183E-01 


0.113346E+01 


-0.421787E-01 


0.305968E-01 


-0.408644E-02 


-0.193771E-01 


0.123304E+01 


-0.270876E-01 


0.185189E-01 


-0.208407E-02 


-0.130319E-01 


0.134340E+01 


-0.166413E-01 


0.107639E-01 


-0.902849E-03 


-0.841563E-02 


0.146641E+01 


-0.978419E-02 


0.600927E-02 


-0.275271E-03 


-0.522831E-02 


0.160437E+01 


-0.557178E-02 


0.322545E-02 


0.208268E-04 


-0.313370E-02 


0.176010E+01 


-0.306684E-02 


0.165657E-02 


0.125739E-03 


-0.181763E-02 


0.193715E+01 


-0.161854E-02 


0.809542E-03 


0.135471E-03 


-0.102257E-02 


0.213998E+01 


-0.824465E-03 


0.375303E-03 


0.110722E-03 


-0.559203E-03 


0.237429E+01 


-0.413458E-03 


0.166244E-03 


0.808655E-04 


-0.297772E-03 


0.264741E+01 


-0.199638E-03 


0.674221E-04 


0.532186E-04 


-0.154652E-03 


0.296895E+01 


-0.923325E-04 


0.240784E-04 


0.325043E-04 


-0.783631E-04 


0.335166E+01 


-0.410173E-04 


0.695348E-05 


0.187892E-04 


-0.386958E-04 


0.381276E+01 


-0.175134E-04 


0.111592E-05 


0.103540E-04 


-0.185709E-04 


0.437596E+01 


-0.741984E-05 


-0.231960E-06 


0.551835E-05 


-0.861492E-05 


0.507462E+01 


-0.295514E-05 


-0.452993E-06 


0.276257E-05 


-0.386774E-05 


0.595690E+01 


-0.110492E-05 


-0.326876E-06 


0.131169E-05 


-0.167290E-05 


0.709434E+01 


-0.387509E-06 


-0.173383E-06 


0.590317E-06 


-0.691226E-06 


0.859696E+01 


-0.126207E-06 


-0.755500E-07 


0.249574E-06 


-0.269811E-06 


0.106411E+02 


-0.375810E-07 


-0.277414E-07 


0.976773E-07 


-0.979770E-07 


0.135237E+02 


-0.999262E-08 


-0.854233E-08 


0.346150E-07 


-0.323958E-07 


0.177768E+02 


-0.228855E-08 


-0.215203E-08 


0.107312E-07 


-0.945155E-08 


0.244325E+02 


-0.428701E-09 


-0.419613E-09 


0.272700E-08 


-0.230841E-08 


0.357133E+02 


-0.637757E-10 


-0.539018E-10 


0.510267E-09 


-0.433381E-09 


0.571587E+02 


-0.579367E-11 


-0.486092E-11 


0.700683E-10 


-0.590489E-10 


0.105975E+03 


-0.256335E-12 


-0.211897E-12 


0.562191E-11 


-0.463011E-11 


0.260189E+03 


-0.281688E-14 


-0.229780E-14 


0.150575E-12 


-0.121213E-12 


0.137016E+04 


-0.682890E-18 


-0.553955E-18 


0.193220E-15 


-0.153367E-15 



where the "step" mass rrix was chosen to depend on the 
number of functions n used in the expansion, so that 

Mo 



1 



(3.50) 



Hence, as n increases, rrix decreases, giving a finer "grid" 
of mass scales. It was found that the precise value of 
Afp was not critical, except that fitting the complicated 
structure of w(2) at large p (to save writing in this section. 



we will sometimes use the notation z[i) to denote the 
generic wave function for model WJC-i, with i = 1,2) 
required a large value of Mq which we chose to be 10 GeV. 
For all other wave functions we chose Afg — 1-2 GeV. The 
leading mass ag determines the asymptotic behavior of 
the wave function at large distances in coordinate space. 



and was therefore chosen to be ao of Eq. (3.441 for the S 
and D states, and rUso from Tables W\ and |V1 for the P 
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TABLE IV. Model WJC-2 momentum space wave functions zi{p), defined in Eq. ( 3.30 1. The momenta are in 
GeV and the wave functions are in units of (GeV)""^". These wave functions are normahzed to 102.3101%. 
The original value of w at the lowest momentum point (-0.743127E-02), which resulted from the imprecise 
cancellation of two very large terms, can be replaced by the bold face number obtained using p^ scaling. 



0.527313E-03 


0.146161E+03 


0.497127E-03 


0.673653E-04 


0.139978E-04 


0.277683E-02 


0.145632E+03 


0.137857E-01 


0.372102E-03 


0.368387E-04 


0.681763E-02 


0.142908E+03 


0.815168E-01 


0.908775E-03 


0.100679E-03 


0.126402E-01 


0.135483E+03 


0.265495E+00 


0.168223E-02 


0.192525E-03 


0.202305E-01 


0.121510E+03 


0.609254E+00 


0.268867E-02 


0.316519E-03 


0.295712E-01 


0.101747E+03 


0.108784E+01 


0.392312E-02 


0.478611E-03 


0.406428E-01 


0.797033E+02 


0.160461E+01 


0.537988E-02 


0.685611E-03 


0.534249E-01 


0.592857E+02 


0.205290E+01 


0.705316E-02 


0.942644E-03 


0.678972E-01 


0.426693E+02 


0.237212E+01 


0.893675E-02 


0.125069E-02 


0.840413E-01 


0.301750E+02 


0.255185E+01 


0.110216E-01 


0.160505E-02 


0.101842E+00 


0.211677E+02 


0.260934E+01 


0.132917E-01 


0.199444E-02 


0.121288E+00 


0.147943E+02 


0.257044E+01 


0.157173E-01 


0.240633E-02 


0.142375E+00 


0.103090E+02 


0.246071E+01 


0.182443E-01 


0.286055E-02 


0.165107E+00 


0.714947E+01 


0.230262E+01 


0.208563E-01 


0.313824E-02 


0.189495E+00 


0.491677E+01 


0.211456E+01 


0.234179E-01 


0.337366E-02 


0.215565E+00 


0.333407E+01 


0.191087E+01 


0.258522E-01 


0.353615E-02 


0.243352E+00 


0.221125E+01 


0.170294E+01 


0.280836E-01 


0.333024E-02 


0.272908E+00 


0.141604E+01 


0.149824E+01 


0.300543E-01 


0.284267E-02 


0.304301E+00 


0.856967E+00 


0.130257E+01 


0.317477E-01 


0.198415E-02 


0.337619E+00 


0.468983E+00 


0.111930E+01 


0.332767E-01 


0.574512E-03 


0.372971E+00 


0.205168E+00 


0.950678E+00 


0.348229E-01 


-0.146236E-02 


0.410488E+00 


0.303889E-01 


0.798954E+00 


0.363889E-01 


-0.389499E-02 


0.450334E+00 


-0.828321E-01 


0.665687E+00 


0.373820E-01 


-0.658983E-02 


0.492699E+00 


-0.154049E+00 


0.551266E+00 


0.364870E-01 


-0.910693E-02 


0.537814E+00 


-0.193834E+00 


0.452057E+00 


0.328528E-01 


-0.113767E-01 


0.585950E+00 


-0.206848E+00 


0.363715E+00 


0.273257E-01 


-0.131172E-01 


0.637426E+00 


-0.198813E+00 


0.284931E+00 


0.211844E-01 


-0.143788E-01 


0.692620E+00 


-0.177269E+00 


0.216739E+00 


0.152062E-01 


-0.150297E-01 


0.751979E+00 


-0.148546E+00 


0.159772E+00 


0.980525E-02 


-0.151078E-01 


0.816027E+00 


-0.118123E+00 


0.113878E+00 


0.537109E-02 


-0.145747E-01 


0.885387E+00 


-0.890513E-01 


0.783167E-01 


0.202329E-02 


-0.134838E-01 


0.960798E+00 


-0.639660E-01 


0.518825E-01 


-0.137523E-03 


-0.119400E-01 


0.104314E+01 


-0.437558E-01 


0.329899E-01 


-0.128366E-02 


-0.101160E-01 


0.113346E+01 


-0.283328E-01 


0.201385E-01 


-0.167665E-02 


-0.815911E-02 


0.123304E+01 


-0.176467E-01 


0.117937E-01 


-0.152450E-02 


-0.629427E-02 


0.134340E+01 


-0.104205E-01 


0.658279E-02 


-0.117834E-02 


-0.462869E-02 


0.146641E+01 


-0.581003E-02 


0.349531E-02 


-0.802520E-03 


-0.324538E-02 


0.160437E+01 


-0.315563E-02 


0.176854E-02 


-0.460812E-03 


-0.218651E-02 


0.176010E+01 


-0.164951E-02 


0.831201E-03 


-0.227186E-03 


-0.141637E-02 


0.193715E+01 


-0.800130E-03 


0.349702E-03 


-0.970097E-04 


-0.879905E-03 


0.213998E+01 


-0.366387E-03 


0.123929E-03 


-0.288285E-04 


-0.527587E-03 


0.237429E+01 


-0.175371E-03 


0.333862E-04 


0.684399E-05 


-0.308052E-03 


0.264741E+01 


-0.775657E-04 


-0.229350E-05 


0.151418E-04 


-0.173993E-03 


0.296895E+01 


-0.303299E-04 


-0.119914E-04 


0.132282E-04 


-0.949967E-04 


0.335166E+01 


-0.104213E-04 


-0.111922E-04 


0.932916E-05 


-0.502335E-04 


0.381276E+01 


-0.309564E-05 


-0.763178E-05 


0.592149E-05 


-0.257338E-04 


0.437596E+01 


-0.144689E-05 


-0.413223E-05 


0.374297E-05 


-0.128285E-04 


0.507462E+01 


-0.488220E-06 


-0.212835E-05 


0.203145E-05 


-0.613856E-05 


0.595690E+01 


-0.107614E-06 


-0.100778E-05 


0.101279E-05 


-0.280394E-05 


0.709434E+01 


-0.588336E-08 


-0.430027E-06 


0.476792E-06 


-0.121466E-05 


0.859696E+01 


0.931326E-08 


-0.163878E-06 


0.211653E-06 


-0.493916E-06 


0.106411E+02 


0.613471E-08 


-0.550345E-07 


0.873897E-07 


-0.185900E-06 


0.135237E+02 


0.246297E-08 


-0.159505E-07 


0.327267E-07 


-0.635345E-07 


0.177768E+02 


0.718045E-09 


-0.386029E-08 


0.106738E-07 


-0.191935E-07 


0.244325E+02 


0.135807E-09 


-0.732832E-09 


0.281290E-08 


-0.493279E-08 


0.357133E+02 


0.340558E-12 


-0.932273E-10 


0.532783E-09 


-0.101591E-08 


0.571587E+02 


-0.140690E-12 


-0.832796E-11 


0.721192E-10 


-0.146914E-09 


0.105975E+03 


-0.105532E-13 


-0.362792E-12 


0.579471E-11 


-0.119175E-10 


0.260189E+03 


-0.149527E-15 


-0.394161E-14 


0.157413E-12 


-0.317980E-12 


0.137016E+04 


-0.319205E-19 


-0.952160E-18 


0.204258E-15 


-0.406280E-15 



TABLE V. The parameters used in the scaling functions 
(02| for model WJC-1. 



TABLE VL The parameters used in the scaling functions 
(021) for model WJC-2. 





N (GeV^) 


R 


m,i (GeV) 


rriso (GeV) 




N (GeV^ 


u 


0.4775 


-0.9856 


1.160 


- 


u 


0.3060 


w 


0.5131 


-0.9942 


0.685 


- 


w 


0.3474 


Vt 


-0.03024 


-0.9772 


1.600 


0.488 


Vt 


0.02700 


Vs 


-0.1038 


-0.9942 


0.840 


0.633 


Vs 


-0.2356 



R msi (GeV) 
-0.9978 1.000 

-0.9805 0.632 

-0.9610 0.700 

-0.9671 1.054 



ruso (GeV) 



0.551 
0.977 
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FIG. 4. (Color online) Ratio of tlie raw wave functions for model WJC-1 (taken from Table [lll| with the probabilities "scaled" to 100% as 
outlined in Table ^ll to the scale functions from Eq. ( |3.42[ l. In each panel the dots are the tabulated wave functions for 60 Gauss points, 
the solid curve the fit, and the dashed curve the fit without the tail wave function of Eq. ||3.47[l. 



states. The mass M^ used in the tail wave functions was 
fixed at 2 GeV for all Zi except u(2), where it was fixed 
at 11 GeV. 

Using these choices, exceUent fits to the momentum 
space wave functions were found. These are shown in 
Figs. [4] and [5] The expansion parameters are given in 
TablesTvnland lVml 

Note that the coefficients of the tail wave functions are 
very small. The tail makes a negligible contribution at 
low momenta, but is very important for the description 
of the high momentum components (for p ^ 3 GeV). 
Remember that the figures show ratios; the actual values 
of the wave functions are quite small above p > 3 GeV, 
as already shown in Fig.lSJ In any case it will be possible 
to study the sensitivity of any observable to the very 
high momentum components simply by setting this last 
coefficient to zero. 



The components z{i) are compared in Fig. p] which 
shows the wave functions for momenta up to 1 GeV, and 
in Fig. [TI which gives an exploded view of the wave func- 
tions from 500 MeV to 1 GeV, a region where all of the 
components might be important, depending on whether 
or not the larger u and w components interfere in a given 
matrix element. 



In general we observe that the P-state components are 
very small, particularly so for WJC-2. However, in the 
region shown in Fig. [7] all of the WJC-1 components are 
larger than the WJC-2 ones, and it is possible that, in 
some observables, P-state components might compensate 
for the differences between the larger S and D-state com- 
ponents. 
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FIG. 5. (Color online) Ratio of the raw wave functions for model WJC-2 (taken from Table [rv| with the probabilities "scaled" to 100% as 
outlined in Table \n\ to the scale functions from Eq. ( |3.42[ l. In each panel the dots are the tabulated wave functions for 60 Gauss points, 
the solid curve the fit, and the dashed curve the fit without the tail wave function of Eq. l |3.47| l. The origin of the many small wiggles in 
u is not clear, but u is very small in this region, as illustrated in Fig. [3] The first Vg point is probably inaccurate because of numerical 
cancellations. 



3. Transformations to coordinate space 



The coordinate space wave functions are constructed 



from the spherical Bessel transforms (3.371. The details 



are given in Appendix C 4 The coordinate space expan- 
sion functions for the terms i < n are 



Q(r) = aA e-^' - e-^' [l + \Z., (l - i?^)] 

1 



G\{r)^AdR,e- 



1 



-Zi 



1 



1 + ^ + \Z, (1 - i?,2) 



G\{T)^B.ARie 



2„-Zi 



where 



-,-Zi 






\zli\-m' 



A,. 



rriir , Zi = 


= Mir , 


^-m; 


mlMf 


B» = 


-A- ^'^ 


{Mf-mlY' 



(3.51) 



(3.52) 
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FIG. 6. (Color online) Momentum space wave functions for both 
models. The largest two wave functions are u (off-scale at small 
momenta) and w (zero at small momenta), with WJC-1 wave func- 
tions solid and WJC-2 dotted. The P-state wave functions are very 
small, with both vt wave functions solid (^'t(l) is negative and vt{2) 
positive) and both Vs wave functions dotted (with {vs{l) negative 
and(Ds(2) very near zero). Note that the S-state wave function has 
a zero around 400 MeV, and that above this momentum the P-state 
wave functions for WJC-1 are no longer completely negligible, as 
illustrated in Fig. [t] 
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FIG. 7. (Color online) Enlarged view of the momentum space 
wave functions for both models. The curves are labeled as in Fig. [6] 
with a dashed zero reference line. The large positive states are ui; 
the large negative ones u, with WJC-1 (solid) and WJC-2 (dotted). 
Both of the vt wave functions are solid, and both Vs wave functions 
are dotted, with WJC-1 negative in both cases and WJC-2 both 
close to zero. 



TABLE VII. The expansion parameters for model WJC-1 wave 
functions. First line is the total number of terms n in the sum 
( |3.48[ ), second is the step mass rux (in MeV), the next lines are 
the coefficients 6f , {i = 1, n— 1}, with the last line the coefficient 
6^, written in the notation 125. (—6) = 125 x 10~^, all in units 
of GeV-^'''^ The "tail" mass M„ = 2 GeV in ah cases. We 
have given the &f coefficients to three decimal places, sufficient 
to reproduce the fits to better than 1% accuracy; if greater 
accuracy is needed use the routine supplied by the authors. 



u 


w 


Vt 


Vs 


17 


16 


12 


12 


75.00 


80.00 


109.09 


109.09 


134.963 


23.813 


-27.120 


-78.763 


52.871 


32.709 


233.155 


688.895 


-217.709 


-111.381 


-998.893 


-2910.298 


1876.699 


844.861 


2679.013 


7690.084 


-11369.449 


-4376.965 


-4894.348 


-13861.596 


49427.176 


16489.297 


6277.553 


17571.825 


-156695.247 


-45122.587 


-5674.481 


-15720.926 


369322.468 


91061.493 


3550.952 


9747.293 


-655367.189 


-136672.748 


-1467.429 


-3994.747 


879178.453 


152411.033 


360.796 


975.086 


-887103.150 


-124633.670 


-40.010 


-107.489 


662786.733 


72617.712 


0.000 


0.000 


-355720.583 


-28549.974 


0.000 


0.000 


129738.381 


6790.953 


0.000 


0.000 


-28805.739 


-738.436 


0.000 


0.000 


2939.797 


0.000 


0.000 


0.000 


-125. (-6) 


-100. (-6) 


498. (-7) 


-394. (-7) 



TABLE VIII. The expansion parameters for each of the model 
WJC-2 wave functions. For an explanation, see the caption to 
Table |Vn] Here the "tail" mass M„ = 11 GeV for u and 2 GeV 
for all other wave functions. 



u 


w 


Vt 


Vs 


30 


18 


21 


18 


344.83 


70.59 


60.00 


70.59 


181.569 


19.228 


33002.947 


-74492.403 


-0.817 


30.600 


-580810.178 


1001357.141 


11.779 


-107.606 


4678729.666 


-6147239.139 


-64.782 


897.969 


-22669425.846 


22510592.546 


247.517 


-5434.481 


72755203.760 


-53617329.732 


-772.915 


25662.551 


-159377601.760 


83750225.006 


1948.950 


-93771.705 


232305154.041 


-77842609.919 


-3813.276 


265664.025 


-193963168.786 


20071103.929 


5397.722 


-582644.240 


15154785.470 


48069105.480 


-4637.491 


986160.259 


160238368.984 


-61037252.628 


533.832 


-1280768.137 


-139567635.494 


8837981.631 


3806.187 


1262888.436 


-52643827.848 


50533251.186 


-3119.875 


-927821.391 


170318699.570 


-65739754.705 


-2058.321 


491712.394 


-69098137.722 


43130588.287 


3773.554 


-177560.697 


-117233319.092 


-16799889.802 


1092.507 


39095.262 


192985378.690 


3717624.361 


-3919.342 


-3959.947 


-138452165.069 


-363261.256 


-736.912 


0.000 


56855870.564 


0.000 


4018.336 


0.000 


-13049893.573 


0.000 


612.635 


0.000 


1310791.818 


0.000 


-4213.343 


0.000 


0.000 


0.000 


-260.753 


0.000 


0.000 


0.000 


4581.603 


0.000 


0.000 


0.000 


-1171.113 


0.000 


0.000 


0.000 


-4559.802 


0.000 


0.000 


0.000 


5566.714 


0.000 


0.000 


0.000 


-2984.961 


0.000 


0.000 


0.000 


814.407 


0.000 


0.000 


0.000 


-92.509 


0.000 


0.000 


0.000 


-118. (-11) 


-172. (-6) 


514.(-7) 


-110. (-6) 
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FIG. 8. (Color online) The coordinate space wave functions for 
both models. The curves are labeled as in Fig. [6] The largest two 
wave functions are the familiar u and w wave functions, with WJC- 
1 solid and WJC-2 dotted. The vt wave functions are both solid, 
with WJC-1 negative and WJC-2 small and positive. The Vs wave 
functions are both dotted, with WJC-1 negative and WJC-2 close 
to zero. 



The last terms with i = n are 
2 

2 

3^ 
2 






MlZlK,{Z^) 



G^{r) = --MlZlK^{Z^), 

lOTT 



(3.53) 



where Kn {z) are the modified Bessel functions of the sec- 
ond kind (see Appendix [C|) . 

At large r the functions G\ [i < n) have the asymptotic 
behavior expected for solutions of the Schrodinger equa- 
tion with orbital angular momentum £. The tail functions 
fall-off like exponentials multiplied by a fractional power 
of r, but, because of their large mass (short range), do 
not contribute to the overall asymptotic behavior of the 
wave functions in coordinate space. At small r the func- 
tions Gl{i < n) have the expected r^'^^ ~ ^i"*"^ behavior. 
However, the tail functions contribute some nonanalytic 
logarithmic behavior at small r, as described in Appendix 
[C| As it turns out, this behavior is too small to be seen 
at the level of one percent. 



4- Wave functions in coordinate space 

The coordinate-space wave functions are shown in 
Figs. [8] and [9J Note that the P-state wave functions van- 
ish beyond 2 fm, and that they are larger for WJC-1, as 
already seen in momentum space. Particularly notable 



FIG. 9. (Color online) Expanded view of the coordinate space 
wave functions for both models at small r. The curves are labeled 
as in Fig. [6] 



is the zero in u(l) near 0.35 fm, leading to a dip at very 
small r. This dip is not an artifact of the fits; it is re- 
lated to the deeper dip in the momentum space u(l) near 
600 MeV, clearly seen in Figs. l6]and[7| A similar (but 
much smaller and not visible in the figures) dip is also 
present in u{2) inside of 0.15 fm. The fact that the dip 
in WJC-1 is larger may be a consequence of its larger 
P-state components. 

Before any physical conclusions can be drawn from the 
differences in the wave functions for WJC-1 and WJC- 
2, observables must be calculated from these models. 
Electromagnetic observables (such as the quadrupole mo- 
ment) will be sensitive to interaction currents, which will 
be different in each case, and might very well compen- 
sate for any difference arising from the wave functions. 
However, these calculations, which will be a subject for 
a future paper, are beyond the scope of this work. 



5. Asymptotic normalization and the D/S ratio 

One observable that can be calculated immediately is 
the ratio of the asymptotic normalization constants for 
the S and D-state wave functions. The results for the 
two models, together with the asymptotic normalization 
As, are summarized in Table [Tx) These quantities can be 



extracted from the expansions (3.51) 



T] = Ad/ As 



Ad = b\"B,Rl 



(3.54) 



Note that rj is more accurately determined than As, and 
that the values of 77 predicted by the two models are 
in agreement with each other (within errors), and also 
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TABLE IX. The asymptotic normalization As and D/S 
ratio rj for the two models. Results are shown for varia- 
tions of n around the central (best) values of 17 and 16 
for WJC-1 and 30 and 18 for WJC-2, with errors based on 
the fluctuations. These may be compared to the recently 
determined r; = 0.0256(4) [T^. 





WJC-1 


WJC-2 


n 


V 


As 


V As 


no -2 


0.02625 


0.86585 


0.02635 0.87471 


no — 1 


0.02623 


0.86837 


0.02639 0.87693 


no 


0.02619 


0.86416 


0.02630 0.87768 


no + 1 


0.02623 


0.86604 


0.02643 0.87829 


no + 2 


0.02622 


0.86436 


0.02614 0.87879 


no (error) 


0.02619(4) 


0.864(2) 


0.0263(1) 0.8777(15) 



agree (to less than 2 standard deviations) with the results 
of Ref. [12]. The agreement is even better with most 
other experimental and theoretical results for r/ quoted 
in Ref. [T^, which are all situated in the range between 
0.0259 and 0.0272. 



IV. CONCLUSIONS 

In this paper we present (Table ll| effective range ex- 
pansions for the ^5*0 and '^Si phase shifts and relativistic 
deuteron wave functions based on the precision models 
WJC-1 and WJC-2 of Ref. I. 

Since the effective range expansions emerge directly 
from the new high precision phase shift analysis described 
in Ref. I, they can be regarded as an up-to-date, preci- 
sion determination of the effective range parameters that 
should constrain any modern theory of the nuclear force, 
such as effective (or chiral effective) field theories. It does 
not matter that they were determined using the CST. 

The covariant deuteron wave functions presented here 
have four components: two that have a nonrelativistic 
analogue, and two (the P-states) of purely relativistic 
origin. Convenient analytic representations of these wave 
functions are presented, with expansion coefficients given 
in Tables VII and |VIII These expansions (and subrou- 
tines that are available from the authors) make it easy 
and convenient to use these wave functions with any cal- 
culation. When these wave functions are used in a co- 
variant theory with one nucleon off-shell they provide a 
precise description of the non-pcrturbative interactions 
that lead to the deuteron bound state. However, even 
in this case a completely consistent description requires 
that we include other contributions, such as final state 
interactions or interaction currents, that are also gener- 
ated by the one-boson-exchange dynamics. That these 
effects cannot be ignored, even in the low-energy or mo- 
mentum limit, is indicated by the size of the derivative 
of the kernel that contributes to the normalization [recall 
Eq. (ISlgl)]. This is about -5% for WJC-1 and -2% for 



energies for some observables. In particular, the calcu- 
lation of any electromagnetic property of the deuteron 
(including the magnetic and quadrupolc moments) that 
does not include interaction currents can be expected ab 
initio to be in error by a few percent. It is therefore quite 
possible that a careful CST calculation of the quadrupole 
moment could resolve the current ^ 5% discrepancy be- 
tween theory and experiment, but this calculation must 
wait until the exchange currents have been accurately 
calculated. 

For this reason it is also unclear how to incorporate 
the wave functions given in this paper into a calculation 
which uses dynamics different from the CST (nonrela- 
tivistic or another form of relativistic dynamics - includ- 
ing light-front). In such cases we are inclined to sug- 
gest that the best procedure is to renormalize the wave 
functions as shown in Eq. (3.40), which is also why our 



WJC-2, and could, of course, be much larger at higher 



numerical wave functions are presented here in this nor- 
malization. 

One observable that does not depend on the normal- 
ization of the wave functions is the asymptotic D/S ra- 
tio. Our results (Table IX) determine this ratio to an 
accuracy about 10 times smaller than the experimental 
accuracy and are in good agreement with measured val- 
ues. Still, it would be interesting to review the theory 
that has gone into the analysis of the experiments that 
determine the D/S ratio and to see if there are any rel- 
ativistic corrections previously overlooked. 
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Appendix A: CONSTRUCTION OF THE KERNEL 

In this appendix we present details, omitted from Ref. I 
[T] , of the construction of the NN kernel directly from the 
sum of Feynman amplitudes that define it. This method 
is completely different from that used in Ref. II ^9^; it 
is simpler and more transparent allowing changes in the 
kernel to be made more easily. 



1. Review of definitions from Ref. I 

In Ref. I the kernels are defined to be 

— 1 f^J PlP2,p'iP'2 (-„ ^ I A ^/J PlP2,P'lP'2 („ \ 

- 2\^dirAiA2,AiAi(P0J + ds»/dirAiA2,-Ai-A^(.P0J 

I c ^r.I P2pl.p'iP'2 ( „ \ I r r^T/^P2Pl,pip2 / „ \\ 

+ °P0 '^dir AiA2,A;A^(."-P0J + <^po<^Sl/dir AiA,,-A; -A^l^-Po)/' 

(Al) 



where 6^^ = (-1)^77(/9t^s)^ (with pT = /OiP2PiP2, A = 
Al — A2 , and / the isospin) is the phase of V under po ~^ 
—po and pi O p2, and Ss = ^p P1P2V (with dp = ±1 
the parity), and in both experssions 77 ~ (—1) is a 
ubiquitous phase. These linear combinations are very 
similar to (but different from) those used referred to as 
Ref II. 

If particle 1 is on-shell, pi — p'l — +. Using the parity 
relation, 



'^dirAiA2,AiA^i=tP0J 



^ TA'/PlP2:P'lP2 ^_l_„ ^ 

PTl/dir-Ai-A2,-A'-A'(±P0J 



(A2) 



the amplitudes (Al) can always be organized so that 
Al = A'l = +1/2, leaving the helicities of particle 2 and 
the phases 6$ and S unconstrained. Hence there are 
2^^ — 16 independent kernels for each J,p2 and P2. Our 
particular choice of independent kernels is denoted Vi and 
defined in Table XI of Ref. I. Their behavior under par- 
ity and interchange symmetry for each combination of 
{/92, p'2} is given in Table XII of Ref. I. 

The evaluation of the matrix elements V"' is simpli- 
fied by the fact that one can write the nucleon spinors 
as direct products of two two-component spinors, one of 
which contains all angle-dependence, while the other de- 
scribes the dependence on momentum and p-spin. Recall 
that the nucleon spinors (from Eq. (El) of Ref. I) are: 



u1ip,X)=Np{pX)(g>xJ0) 



ug(p,A)=iV,(pA)$5x_.(^) 



u(P,^) P = + 
w(-p,A) p=- 

|m(-P, -A) p^ + 
|i;(p,-A) p=~ 



where the p-space spinors are 

cosh ^( \ 

2AsinhiC/ 

-2AsinhiC 

cosh ^C 




(A4) 



with p = IpI and tanh^ = p/Ep. Note that, at large p, 

cosh ^( 
sinh IC 




yJ1m(Ep + 771 



(A5) 



For momenta limited to the xz plane, the spin 1/2 spinors 
are 



X,„(0) = i?,(0) 



X^.I2%^RyiO)\ \ 



sin \B ) 
- sin \6 

cos ifl 



(A6) 



2. Extracting the angular integrals 

Each term in the Feynman meson-exchange operator 
can be written in the generic form 

Vb(0,po) = ^5^s,2?b(0,po) 

ij 

x[0^(j>^)®S^]\[0\p^)®S=]\, (A7) 

where O* are operators in the p-spin space, S^ opera- 
tors in the spin-space, gi and Sj are constants, and the 
sums are over all operators needed to describe the meson 
exchange interaction. Here I?^ is the scalar part of the 
propagator 



T^b{0,Po) 



fbi^b,q) 



ml + \q{9,poW 



(A8) 



with q{0,p(j) the four-momentum transferred by the ex- 
changed meson and /;, the meson form factor. The oper- 
ator 0{pq) acts in the 2x2 p-spin space, and S operates 
in the 2x2 spin space, so that O (E) S is a, AxA matrix. 
Using this notation, the matrix elements of V reduce to 
(suppressing the index b for simplicity) 



yPlP2,PiP2 ^ 

dir Al A2 : A^ AJ, 



-1-^ ^ _^OPlPl,P2P2 



(A3) 



(AiA2|A'iA^)I?(0,±Po). (A9) 
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where 



of the form 



O 



P1P1.P2P2 



i 



(AlO) 



with N = N^T'^ the Dirac conjugation in p-spin space, 
and use is made of the fact that the partial wave amph- 
tude may be calculated by aligning the initial momentum 
p' in the +z direction. Note that the sign of po in the 
operators 0{±Pq) is captured in the last subscript of the 
matrix elements O, that the superscripts for O group 
together the p-spin indices for particle 1, followed by the 
p-spin indices for particle 2, and that the matrix elements 
of S do not depend on the p-spin of the states. We keep 
Ai unspecified in these formulae even though only the 
case Ai = + needs to be considered. The linear combi- 

-, but the individual terms in 



nations ( Al I restrict X'^ 

the sum use both signs of X'^ , so it cannot be restricted 

in (AlO). Recalling that the partial wave projections are 



^di^A^r.lfA^bo) =2n f\in0d9di,,{0) 



^virXtx'SP^P'-^n (All) 

we see that each of the matrix elements (AiA2|A']^A2) will 
be multiplied by the rotation function dy^{d), where 
A = Al — A2 and A' = A'^ — Aj. It turns out that the 
/9-spin matrix elements generate an additional factor of 
z — COS0, so matrix elements with an additional factor 
of z are needed; these are handled by multiplying each 
of the matrix elements by z% where i = or 1. The 16 
integrals we require are defined in Table pC} with explicit 
expressions for these 16 matrix elements in terms of 5 



independent integrals summarized in Table XI The 5 in- 
dependent integrals over Legendre polynomials and their 
derivatives are given in Table |XII[ with convenient iden- 
tities showing how to reduce 5 other integrals to these 5 
given in Table pOnl 

The 16 generic integrals (for each value of pi) defined 
in Tables XI and XII of Ref . I can then be written as prod- 
ucts of the 16 angular-dependent terms given in Table [X] 
and Dirac matrix elements (defined below) as follows: 



,,P2P2 _ 



,,P2P2 



,,P2P2 



P2P2 
'' 9 



P2P2 



,,P2P2 _ 



„P2P2 



E 

i 

E 

i 

E 

i 

E 

i 

E 

i 

E 

i 

E 

i 

E 



^/^Pl + -P2p'2,i Ai _ ^/-iPl+,p2p'2-i Ai 



o 



Pl + ;P2p'2,i Ai /oPl + :P2P2,'' Ai 



/■/^Pl + :P2P2:i Ai 



AlPl + .P2P2.« Ai 



-,pi + .p2P2,i Ai 



/OPl + ,P2P2:* Ai 



^/^Pl + ,P2P'2-i Ai 



Pl + .p2P2-i Ai 



/-/1P1T,(>2P2>' 41 , /■qtJl-r.f2IJ2'^ At 



/r^Pl + -P2P'2,i Ai , /■/^Pl+,P2P'2,i Ai 



/y-,Pl + ;P2p'2,i Ai 



^/^Pl + :P2p'2,i Ai 



iy-)Pl + ,P2p'2,i Ai 



C^Pl + ,P2P2A Ai 



± 



± 



± 



± 



± 



± 



/-/lP2 + ,PlP2 Ai 



/-/-)P2 + ,PlP2 Ai /^P2^,P1P2 Al 



_ qP2+,P1P2 Ai 



-)P2 + ,PlP2 Ai 



//^P2 + ,PlP2 Ai 



/■/^P2 + ,PlP'2 Ai 



//-lP2 + :PlP2 Ai 



^P2 + ,PlP2 Ai 



/-nP2 + ,PlP2 Ai , //^P2+,PlP2 Ai 



//-lP2 + ,PlP2 Ai I /■/^P2 + ,P1P'2 Ai 

'-^-l— , + - ,-^-1— ,-|— +'-^ + -,- + ,-^-1- 



/nP2 + ,PlP2 Ai 



//lP2 + ,PlP2 Ai 



//-)P2 + ,PlP2 Ai I ^/-,P2 + ,P1P'2 Ai 



(A12) 



r 



where the -I- and — integrals include the propagator 

= 7r/ sm0d9V{0,±po) (A13) 

Jo 



and, in each case, the upper sign goes with the upper 
index (of the two lower indices) and the lower sign with 
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TABLE X. The 16 matrix elements A 
is written as A = ± and z — cos t 



AiA2,A'^A2 



= z' (+A2lA'iA'2) di,^ for Ai = + and i = or 1. 



(+ 
(+ 



+ 
+ 



^> rfoo 
h) did 
-) d(o 
-)d-L 



z\\{l + z)di^ 
z" \ sin Q doi 
-2:'|sin6ldfo 

Z'|(l + 2)dfl 



A\ 
A\ 
A\ 



+ 



+ - 



+ 



Here A = ±| 



-) 4. 

-) dill 
+) d-i, 
+) d'-.. 



J 1 



~! 1 



(1 - z)di 



■00 



z' -^ sin ^ dpi 
z^ I sin 61 d'Li q 
z^l{l-z)dU^, 



and og from Table 
needed only for the 



XIII 



Note that d:iio — ~diQ and dfji = — dio- Here a^ = ai{9). 
and hence only for vector mesons. The matrix elements Aij,k 
the code 1 = ++, 2 = +-, 3 = — h, 4 = and fc = 1 for A° and fc = 2 for A^ . 



TABLE XL Matri x elem ents of Table [X] expressed in terms of the angular functions from Table [XlT] and ag, 07, ag 
— A „_ j-^„™ 't^„ui„IvtttI m„+„ +!,„+ jJ _ j^ „„j j~r _ j^ u Jia\ 'pj^g integrals 04 and 05 are 

have i,j — {1,4} according to 



Al..^. 



I-+ 



= ^11,1 
= ^11,2 
= ^21,1 

~ .421,2 
= .4l2,l 
= .4l2,2 
= .422,1 
= .422,2 



|(l + 2)d^0 = 
\z{l + z)dio 

I sin 9 d'oi — - 
dii = 
' I sin 6 d'lo - 
' ^ z sin 6 d'ld 



|(ai +02) 
= |(a2 + (14) 

1 , 



jag 



jflg 



= |(l + 2)dfi = i (as + ai) 
= |z(l + 2:)dfi = |(a9 + a2) 



AU 


__ 


= .4l4,l = 


-|(1 - 2:)d^o = -|(ai - 02) 


AU 





= .4l4,2 = 


-|2:(1 - z)dio = -|(a2 - a4) 


A^ 





= A24,l = 


1 sin ^ did ~ ^'^<^7 


^i- 


__ 


= .424,2 = 


^zsinSdid = — fag 


A% 


— h 


= .4l3,l = 


1 sin0d£.i_o = — |a7 


A+ 


— h 


= .4l3,2 = 


lasinSdiio = —-^0,8 


A'+- 


— h 


= .423,1 = 


|(l-2)d:ii,i = |(a6-ai) 


Ai- 


— h 


= .423,2 = 


|2:(l-z)d£i,i = 1(09-02) 



TABLE XII. The five independent angular 
integrals, a'' (6), needed for this calculation. 
Four more frequently used linear combina- 
tions of these, denoted Oj, with j = {6,9}, 
are given in Table [XlII] 



ai{e) 

ai{e) 
ai{e) 



zPj{z) 
Pj-i{z) 
z^Pj{z) 
zPj-i{z) 



the lower index (of the two lower indices), so that, for 
example, vi is symmetric under po ~^ ~Po and U5 is an- 
tisymmetric. The 16 p-spin matrix elements that enter 
into these expressions, all independent of 9, result from 
the fact that the general /9-spin matrix elements O given 
in Eq. (AlO) are linear in z, and can therefore be sepa- 



rated into two terms using the notation 



O' 



PlPl,P2P2 

Ai A2,A^ A2,it 



o 



Plp'l,P2P'2,0 

Ai A2,A^ A2,i 



-0^:t^^- (A14) 



The O^'s arise only from the exchange of vector mesons; 
all of the Cs arising from the exchange of each meson 
will be computed below. 



TABLE XIII. Useful identities involving the 
d functions. Here Pj — Pj{z), ai = af{6), 
and the new combinations, a6,a7,ag, and ag 
are defined. 



I [(1 + z)did - (1 - z)dU,{\ = dio = Pj 



ffli 



\ [{l+z)d-h + {l-z)d-U^{\ =06 

= -j^ [JzPj + P.j_.i] = -j^ [Ja2 + aa] 

\z [(1 + z)did + (1 - z)d:L^^^ = ag 

= T+i [Jz^Pj + zPj-i] = j^ [Jai + as] 

sin^d^'o = ^7^ [zPj - Pj-i 

J [02 - as] = 07 



z sin 9d' 



■J 




04 - as = ag 



3. Particle momenta 

In some applications there is the possibility that both 
particles in the final state could be off-shell (the possibil- 
ity that both particles in the initial state are off-shell is 
not considered at this time) . To allow for this possibility, 
we will write the final state momenta in the cm system 
of the two-nucleon system in the general form 



Pi 

P2 



{~xoEp + {l + xo)lW,-p} 



with 



Pa = xoiEp - ^W), 



(A15) 



(A16) 



where xq is a dimensionless number varying between —00 
and 00, and W is the total two-body energy. Note that 
when a;o = 1, particle 1 is on-shell (with energy Ep), while 
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when xq = —1, particle 2 is on-shell. Hence changing 
the sign of xq is a convenient way to interchange the 
energies of particles 1 and 2 in the final state, and we 
may construct the ±po combinations of Eq. ( Al I merely 



by changing the sign of xq. (A similar variable, yoi could 
be used for the initial state.) 



4. Single particle p-spin matrix elements 



We now turn our attention to the form of the p-spin 

matrix elements C'^^^^'vv'*- -^^ ^^^ ^^^ ^^^ simplest 
cases, these are best calculated by first constructing ma- 
trix elements on each nucleon line and then multiplying 
these together to get the total O. The one-nucleon matrix 
elements are evaluated numerically by matrix multiplica- 
tion, best described separately for each type of meson ex- 
change. The total p-spin matrix elements for each meson 
will be summarized in the next subsection. 



a. Scalar mesons 

The calculation of the on-shell scalar exchange is very 
straightforward, with 

-)PlPi,P2P2.0,S 



O 



Al A2,A-|^ A2 ,i 



Al A2,Aj^ A2,i 



Np,{pX2)lNp.Jp'X'^) 



(A17) 



where 1 is the unit matrix in 2x2 space, and the addi- 
tional superscript s labels these Cs as the contributions 
from scalar mesons. It is convenient to calculate the gen- 
eral matrix element 



Np^{pX,)lNp,{p'X'^ 



(A18) 



and then construct all of the products ( |A17 ) , avoiding as 
much algebra as possible. This also allows the off-shell 
matrix elements to be calculated without any extra work. 
In what follows we will generalize the notation of (AI8I: 



O^ 



NpSpX,)ONp.^{p'X',) 



(A19) 



where O is any 2x2 operator in the p-spin space. Explic- 
itly, the Dirac matrices are written as a direct product 
of two 2x2 matrices (where the first matrix operates in 
the 2x2 p-spin space and the second in the spin space): 



7^ ^ -r^ (g) 1 



7 = T (g) 1 

7» = iT2(g)cr^ (A20) 

Hence, the p-spin matrices in the 2x2 spin independent 
part of the Dirac space are the familiar Pauli matrices. 

We frequently encounter off-shell couplings, which in 
every case give factors of 

(m- f) uP{-p, A) = (to -h 7 • p - 7%) uf (p, A) 

= (pi?p-po)7V(p,A) (A21) 



where we have used the Dirac equation for the (always) 
on-shell spinor u''(p,A) and recalled that M^(p,A) = 
w(— p. A) so that its Dirac equation is 

{m+ ^)v{-p, A) = (m + j"Ep + 7 • p)w(-p. A) = 0. 

This use of the Dirac equations allows us to replace the 
angular dependent 7 • p term by an angular independent 
factor. (Similar steps work off-shell couplings in the final 
state.) Off-shell couplings therefore do not involve any 
new angular integrals, but they do require evaluation of 
a new 2x2 Dirac matrix element. Since the most general 
scalar Feynman operator is 



^s{p,p') = gsl-i^s 



2to 



2to 



(A22) 



the spin independent part of the matrix element in the 
general scalar case is then 

A,,{p,,p'i)=uP'{p,X,)As{p,,p'^uP'^{p',X'^ 

= 9s 1* - ^-rf [p,Ep - pm + p[Ep, - p[^] 

^R';:j:,j±po,±p'o), (A23) 

where, in the last line, we define the first of the one- 
particle p-spin matrix elements needed in this calcula- 
tion. The definition of Rs restores explicit reference to 
the helicities A^ and A,^ implicitly contained in the ma- 
trix elements 1^ and r?, and replaces pio by ±po, with 
Pio = W/2 + po and p2o = W/2 — po (so that sign of 
Po is positive for particle 1 and negative for particle 2). 
We emphasize that the spin dependent angular part is 
unchanged by the off-shell couplings. 

We encounter the off shell factors numerous times, so 
it is convenient to denote 



AEi = piEp -pio 

AEl = p'.Ep, - p[o. 
However, note that 

AE1AE2 = AE[AE'2 = 
because one of the two particles is always in shell. 

b. Pseudoscalar mesons 



(A24) 
(A25) 



The pseudoscalar Feynm 


in operator is 


%{P,P') ^9pl^ 


-5p(l-Ap)| 


m— ^ 


f +f 


2to 


2to 



(A26) 



Using the same arguments, the pseudoscalar matrix ele- 
ments become 



KiiP^^P'i) ^9p\ tI 7, — — i'^f 

[^ 2m 



AE, - AE' 



(A27) 
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TABLE XIV. Matrix elements (p, A|0„|/5', A') = Np{pX)OnNp,{p' X') = -DJ^, „ for various operators On, 
where n = {1,4} with the correspond ence r^ — >■ 1, ir^ — >■ 2, r^ — >■ 3, and 1 — >■ 4. The two-component 
p-spin spinors Np were defined in Eq. (A4l. In this table c = cosh |(^, c' — cosh|("', s — 2Asinh|(", and 



s' = 2 A' sinh jC', where ^ was defined below Eq. ( A4| 



(+,A111+,A') = - 
( + ,A|r^|+,A') = 
U,\\t^\+,X'} = 
{ + ,X\ir^\+,X') = 



{-,A|1|-,A') =cc' -ss' 
-,X\r^\-,X') ^cc' + ss' 
-,A|r\A|-,A') ^cs' -sc' 
-{-,X\ir^\-,X') =cs' + sc' 



(+,A|1|-,A') = (-,A|1| + ,A') = -{cs' + sc') 

{ + ,X\r-'\-,X') = - (-,X\r'\+,X') ^ -{cs' - sc') 

+ ,A|rV>^') = -{-,Mt^\+,^') ^cc +ss' 
+ , Alir^l-, A') = (-, A|ir^ A| + , A') = cc' - ss' 



where we introduced the new matrix element Rp. This 
requires evaluation of the matrix elements of r^ and «t^, 
but the spin dependent matrix elements are the same 
as in the scalar case. For the evaluation of the matrix 



elements, see Table XIV 



c. Vector mesons 

The vector matrix elements introduce new spin depen- 
dent factors, and the terms O^ linear in z. 

The first step in the reduction of the vector-meson ex- 
change terms is to reduce the Pauli interaction term using 
the generalized Gordon decomposition: 



2m 



= r 



(p + pT 

2m 

TO— ^ 

2to 



^(^ - ryt^ 



m— ^' 
2m 



(A28) 



Combining this result with the general definition of the 
vector Feynman operator 



K{p,p')^g.W + v, 



m— ^ 
2m 



7^ + 7'' 



2to 



+ ^i.-(.-p'). 



(A29) 



gives 



A^'(p,p') = avW[\ + K.) - ^kp + p'Y + (^. - ^.) 



m— ^ 
2m 



7^' + Y 



2m 



^,L.a + ».)_|i,,+,y + l^ 



b.Ef^^' + r^i'rflS.E' 



(A30) 



where the last line anticipates later use of the Dirac equation to reduce the operators m— ^. 
These two vertex operators are contracted with the spin-one meson propagator, which gives 

A^i(pi,p'i)A:2b2,P^)A^.=A^(pi,p;)AS2(p2,p'2)-A^„i(pi,p;)A^,2b2,Pi) 

"^ ^{At'i(pi,pi)(pi -p!\)i,K2{Pi^A){P2-p!2)y\ 



(A31) 



where the factor 77 = 1 is included solely to keep track of the effect of the q^(f jn?^ in the meson propagator (note 
that the sign follows from —q^q" — (pi —p'iY{p2 —p'2Y)- To simplify, note that 



A",fe,p-)=5JTf(l-fK„) + l« 



|^(Ai?,-)-A£;:)-|^(p,i?, + p^i?,, 



i?^f^V.i(±Po,±Po), (A32) 



r 



where we introduced the new matrix elements i?„i and Rvi- 
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The three-vector part introduces the spin-dependent 
operators cr' and angle-dependent terms from the factor 
(p -f p')\ These are first separated into two terms, 



A',,fe,pD = A, xd^' 






-P'. 



'V 



(A33) 



where we remind the reader that the subscript i labels the 
particle number (1 or 2), and Ai is the spin-independent 
coefficient of the spin-dependent operator a^ , and the 
spin-dependent part of the second term is the identity, 
suppressed in the expression. The coefficient Ai is 



A,^g^^r^{l + K,)+rl^^^^^-^{^E,-^E[) 



2m 



<v,.3(±Po,±Po)> 



(A34) 



which defines the matrix element Ry^. 

In calculating the operator A^;^(pi,p'^)A:^ 2(^2,^2) ^6 
encounter the cross terms 



(Tl ■ (P2+P2) 



-CTl • (pi +p[) 

-2(Aip + A'ip'), 



(A35) 



where the evaluation has been carried out in the center 
of mass system and use was made of the fact that the 
helicity eigenstates satisfy the eigenvalue conditions 



cr-p|A) =2Ap|A) . 



(A36) 



(For inconnng particle 2, for example, 2) is a state with 
spin down in the z direction, and pj is in the —2 direc- 
tion, so that <T • P2 |i) = —pa^ |i) = pj^).) Finally, 
note that 

(Pl + Pi) • (P2 + P2) = -(Pl + Pi)' 

= -(p2+p'2 + 2pp'z), (A37) 

where we use p^, p^, and p to denote the square of 
the three-vector, the four-vector, and the magnitude of 



the three- vector, respectively. Combining (A33) - (A37) 
gives the following: 

Ki{Pi,Pi)K2{P2,p'2) = ^1^2 (AiA2|cri • cr2|A'iA^) 
+ I ^^^ Uil2(Aip + A'lp') + A2li{X2P + \y) 

771 L 

-(^)\il2(p2 + p'2 + 2pp'z)|(AiA2|A;A^), 

(A38) 

where, for clarity, we have included the term (AiA2|A']^A2) 
that multiplies the terms in the {}; this term is not part 
of the p-spin matrix element and will be removed once 
the effect of the operator cri ■ cr^ has been expressed using 
the relation 



(AiA2|cri • cr2|A'iA;,) ^ r (AiA2|A'iA 



2/ ; 



(A39) 



where T will be determined shortly. The extra z ~ cos B 
dependence in the last term of ( A38 ) is the origin of the 



special linear z dependent terms described in Eq. (A14) 



above. In our model, only vector-meson exchanges con- 
tribute such terms. 

To determine T correctly note that the helicity states 
I -I — ) or I — h) correspond to the spin states |tt) or jH), 
respectively, and hence are spin triplet states with T = 1. 
The same argument works for both initial or final states. 
This leaves only the combinations with equal helicities 
in both the initial and final states: Ai = A2 = ± and 



-^i = ''^2 = i- Looking at Eq. (A12), the only potentials 



affected are wi,W5,W9, and V13, and using the results 

cri-cr2|++) = -|++)+2| ) 

Ti-<T2|— > = -|— )+2|++) (A40) 

the T term in these kernels leads to the decomposition 



,,P2P2 



^/^ + + ,P2P2.«-T- /. i _|_ (^ 



+ +,P2P2,*T- \ i 






hs{oti:'2^t T ot^'^^n ( A^. 



'++,+^ 



h+, ,+ 



+ ,++ 



A' 



A\ 



+, — ) ' 
(A41) 



where S = —3. If the T operator is absent, the same 
expression holds, but with S" = 1, so this decomposi- 
tion can be used for both cases. With this relation it is 
not necessary to explicitly calculate the matrix elements 

(AiA2|<Ti.^2|AiA'2). 



d. Axial vector mesons 

Only the simplest possible coupling was used for the 
axial vector mesons, namely 

K^(P^,P^ - 9alH , (A42) 

and this was contracted with the g^i, tensor, giving 



A^i(pi,p'i)A:^2(P2,p'2)5m. = A°iA02 (AiA2|A'iA'2 

— C1C2 (AiA2|<Ti • <T2|AjA2) , 



(A43) 



with 



Al,^~tgar!^Rl]''l^,{±Po,±p'o) 

C,^-gaT!^Ri]''l^^{±Po,±p'o), (A44) 

which defines the last two p-spin matrix elements, Rai 
and Ra2- 

5. Two-particle p-spin matrix elements 

The two-particle p-spin matrix elements are the ones 
defined in Eq. (A14|. They are constructe d from the he- 
licity matrix elements C^y ^ given in Table . 



XIV 



and the 



one-nucleon matrix elements Ri, defined in Eqs. (A23l, 



(A27), (A32I, (A34), and (A44) above. In constructing 
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the matrix elements, note that po — > — po corresponds to 
the mterchange of pio O p2o, and similarily for Pq. The 
resuhs, with all indices included, are: 



Scalar, from Eq. (A23) 



O' 



PlP'l,P2P'2,0,S 

Ai A2,A^ A2,i 



-R 



i7a^po,p'o)RZ%Mpo^-po) 






(A45) 



Pseudocalar, from Eq. (A27) 



O' 



AiA2,AiA^,± - ^AiA;,p(.±PO,PoJ^A2A;,,p(=FPO,-Po) 



/-r)PlPi.P2P2>l,P _ Q 

Ai A2 ,Aj^ A2,ib 



(A46) 



Vector, from Eqs. (|A32|, ( |A34[ ), and ( |A38[ ): 

C^Plp'l:P2P2:0,l' _ „pip'i f + nn „MK-''2P2 ('zr„„ _„M 

*-^AiA2,AiA^,± - -"AiAi,t,ll='=PO,PoJ^A2Ai,ul'>TPO, PoJ 



77 



i?' 



r/AU2(±P0,P;))i?rA!,.2(TP0,-p^) 



-rKX^,{±Po,p',)RZt^3iTPo,^p'o) 

+ {X,p + Kp')RZt.,iTPo, -p'o)DZX,_ 

(A47) 



AiA2,AjA2,± -f^^ \^ 2m 



7-)PlPl nP2p2 

^AiA;,4^A2A^,4- 



In these expressions we used the fact that D^ does not 
depend on po ■ 



As a first step, define artificial couplings designed to 
reproduce the ratio G'2/Gi ~ Koff = i^p = 1.7930 and 
gl^ = Gil Up = G\IGi = -0.00779. The upshot of this 
is to use the proton anomalous moment for k and a ficti- 
tious charge of g\^ = —0.00779 so that g^ff'^p = ^1 ^"^^ 
(7gfj Kp = G2 are reproduced. Then, we will factor out the 
coupling (7cff from the i?'s so that it occurs as an overall 
multiplicative factor with the correct (negative) sign. 

Next, note that all terms in the photon exchange must 
involve at least one factor of k, so that B ^ R^^ does 
not contribute. The i?„i and _Rt,3 squared terms must 
have the g\^ terms subtracted. Noting that i^ = for the 
photon, and redefining the photon R 's so that g^ — 1 and 
Kv = Kp, this gives 

Photon: 



ipipi 



^AiA'i,?^ 
2 f r>PlP 



Ai A^ ,3 A2 A2 ,3 
' 5cff \-"^AiAi,73 



i±PoA)RZ''l.3 



{TPo,-Po) 



— D''^''^ JjP-^P-2 



Go 



-D 



piPi 



.D 



P2P2 



Gir 

m 



Aff-f^ AiA^,4 A2A2,4 



{X^P + X[p')RZ''l,3i±Po>Po)DZ%A 
+ {X,P + X',p')RZ%,,{TPo,-p'o)Kka 



o 



PlPi.P2P2:l:7 _ /-I PP T-fPlp'l T^P2P'2 

AiA2,AiA;j,± - ^2 2^2-^AiA;,4-^A2A^,4' 



(A50) 



where the redefined i?'s are [recalling Eqs. (A32) for Ryi 
and ( A34 1 for Ry^ with v ~ and g^ = 1] 



Axial Vector, from Eqs. (A43l and (A44) 

Plp'i:P2P2iO,a 



o 



o 



Ai A2,A-]^ A2 ,i 

pip'i,p2p2A,a 



d 



{^zi. 



£)P2P2 

2 A2Ar,.2 



■-f£)P^Pl JjP^P 
Ai At ,3 A2 A 



0. 



(A48) 



Ai A2,Aj^ A2,i 

Once again we note that the D's do not depend on po 



^rA',7i(^o,2/o) = (i+A^p)^r:A',3 



^DZ,,ApEp + p'Ep,) 



RZ',^3i^o,yo)^{l + ^p)DZr 



"■p jjpp 
2m ^^'^1 



AEp{x„)-AE'{yo) . (A51) 



Treatment of the photon 

For the photon we assume that Fi„ = and that all 
other form factors are equal: -P^n = Fip = F2p = Fo- 
The photon is an isovector-vector exchange, obtained us- 
ing the following couplings (with a ~ 0.007297): 

(7 — > no e term 
g^K -^ e^K^ -> -a(1.913) = -0.01396 = Gi 
-a(1.913)(1.793) 



2 2 
g K. 



i-' i^lt f^Ti 



-0.02503 = G, 



(A49) 



These couplings require a special construction and cannot 
easily be incorporated into the general vector formulae 
above. 



6. Isospin factor 

For the exchange of isovector mesons, one must evalu- 
ate the factor Ti • T2 . This operator preserves the isospin 
symmetry of the state; therefore its value is the same for 
each of the 16 kernels (or 4 when J = 0) that make up 
the coupled array describing the singlet (S), triplet (T), 
or coupled (C) channels: namely it must equal unity on 
isovector np states and —3 on isoscalar np states. To de- 
termine the correct value look at the leading amplitudes 
(those which are symmetric under the change in sign of 
Po, and hence survive when both particles are on shell), 
and note that, for these amplitudes, the Pauli principle 
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requires that L + S + I he odd. For singlet and triplet 
states, J = L while for coupled states J ~ L±l. Hence 
we have the identification: 



singlet & coupled 



done) and the first term arrises from vr*^ exchange and 
the second from tt* exchange. Noting that 



triplet : 

These are summarized by: 

singlet & coupled : ti • T2 = 
triplet : Ti • T2 = 











t!tI \pn} = 


= - \pn) 


J even 


-> / = 


= 1 




t!tI \np) = 


= - \np) 


J odd 


^ / = 


= 




2(^1+^2"+ ^1+^2") H> = 


. 2 \np) 


J even 
J odd 


^ / = 

-^ 1 = 


= 
= 1. 


(A52) 


2{T+T^+T+T^)\np) = 

we see that 


-- 2 \pn) 



(A55) 



-l + 2(-l)^ 
-l-2(-l)'^ (A53) 



To determine the correct isospin factor for each t'f^''^ 
the channel to which it contributes must first be deter- 
mined by examining the matrices Eqs. (E43)-(E48) of 
Ref. I. Inspection shows that the kernels wi,W3,W9, and 
I'll do not contribute to the triplet channels, while ker- 
nels i'5,W7,wi3, and U15 contribute only to triplet chan- 
nels. The rest contribute to either, depending on their 
/9-spin. In particular, contributions to the triplet chan- 
nels come from 

• the ++ sector of vq and v^, 

• the H — sector of viq and t;i2, 

• the — h sector of V2 and W4, and 

• the sector of vn and vie, 

with the other sectors of these kernels contributing to 
either singlet or coupled. 



7. Charge symmetry breaking induced by pion 
exchange 

In order to include important violations of charge sym- 
metry and to better describe the important pion ex- 
change, both models treat the charged tt* and the neu- 
tral 7r° as independent mesons, with their masses fixed 
to their experimental values. Model WJC-1 allows the 
couplings of the pions to be varied independently while 
model WJC-2 assumes that the couplings (written in 
terms of g and not /) of both are the same. 

In addition to the couplings and masses, the correct 
isospin factor for the exchange of each pion must be 
worked out. There are two ways to do this. The op- 
erator Ti • T2 can be decomposed into 



Tl ■T2 



rM 



'■{r^ 



'1 '2 
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where r = (T^±r^)/2 are the raising and lowering oper- 
ators (normalized to 1 and not y/2 as is more commonly 



Tl -1-21^0 ^ -1 



Tl -1-21^* -^ 



2 7 = 1 states 
—2 7 = states . 



(A56) 



These factors can also be obtained by using the fact that 
isospin invariance gives couplings for {p — > pn'^ , n — > 
mT^,p — >■ mr~^,n — >■ p7 r~} i n the ratios {1, —1, \/2, \/2}. 

The impHcation of ( A56 ) for tt^ couplings is to replace 
the isospin factors from ( A53 ) by 



Tl ■T2\^± -^ 



2(-l)-'' Sand estates 
-2(-l)'^ T states 



(A57) 



Appendix B: EVALUATION OF ANGULAR 
INTEGRALS 



1. Angular integrals 

The angular integrals needed for this problem are de- 
fined to be 

(Bl) 

where Vf, is the dressed propagator (including form fac- 
tor) for meson b, defined in Eq. (A8), and the four- 
momentum transfer for the direct (+) and alternating 
(— ) terms is 



q^ixo) = [xoiEp - iM^) - yo{Ep, - ^W) 
-p2 - p''^ + 2pp'z, 



(B2) 



where in this paper the initial state always has particle 
1 on-shell, so that yo = 1. These integrals are assembled 
from the matrix elements of Table |X] and the functions 
of Table pOTTl leading to the identities of Table |Xl| All 
of the integrals A are recalculated for each J,p, xo,p', 
and yo, so these arguments are not included in the list of 
subscripts denoting the dependences of A. 



2. Mapping for the direct terms 

For the direct term the meson propagator is strongly 
peaked in the forward direction (z = 1). For the recent 
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TABLE XV. Convergence of the diagonal element of the direct part of the kernel v^^ at the point p 
use no mapping; those labeled z{y) use the mapping function ( |B5| . Here the 
factor, all of these matrix elements are reduced by a factor of 7xl0~^^. (All 

J = J = l J = 2 J = 



■ P = Pmax- Columns labeled y 
nucleon form factor has been set to 1. With the correct form 
numbers in this table are multiplied by 10~^) 

3 J=4 J=5 J= 



6 



n y z{y) 


y z{y) 


y z(y) 


y z{y) 


y z{y) 


y z 


(y) 


y z{y) 


20 0.14321 -3.3764 


-0.38372 -24.288 


0.14001 -3.3825 


-0.36991 -24.233 


0.13277 -3.3967 


-0.34585 -24.134 


0.12190 -3.4184 


30 0.62029 -3.3765 


-2.0149 -24.288 


0.61385 -3.3827 


-1.9816 -24.232 


0.59907 -3.3968 


-1.9225 -24.133 


0.57642 -3.4186 


40 1.4365 -3.3765 


-5.2403 -24.288 


1.4277 -3.3827 


-5.1898 -24.232 


1.4074 -3.3968 


-5.0999 -24.133 


1.3762 -3.4186 


60 2.5656 -3.3765 


-12.496 -24.288 


2.5569 -3.3827 


-12.435 -24.232 


2.5368 -3.3968 


-12.326 -24.133 


2.5059 -3.4186 


90 1.4055 -3.3765 


-18.312 -24.288 


1.3992 -3.3827 


-18.254 -24.232 


1.3847 -3.3968 


-18.152 -24.133 


1.3624 -3.4186 


TABLE XVI. Convergence of the diagonal element of the exchange part of the kernel v^ at the point 


p=p' =Pma,L- 


Columns labeled 


y use no mapping; those labeled z{y) use the mapping function 


(|B8|l. Here the nucleon form factor has been set to 1. 


With the correct 


form factor, all of these matrix elements are 


reduced by a factor of 7x10 ^'^ . 










J = 


J = l 


J = 2 


J = 3 


J = 4 


J=5 




J = 6 


n y z{y) 


y z{y) 


y z{y) 


y z{y) 


y z{y) 


y z{y) 


y 


4y) 


20 -25.184 -74.786 


-5.1601 -14.525 


-16.551 -49.079 


-2.1999 -6.2065 - 


2.1638 -6.2802 0.98501 2.7617 


8.5080 25.368 


30 -37.148 -81.841 


-6.0514 -11.903 


-24.393 -53.713 


-2.5703 -5.0887 - 


3.1487 -6.8793 1.1690 2.2584 


12.580 27.757 


40 -43.422 -80.907 


-8.1636 -13.401 


-28.504 -53.097 


-3.4750 -5.7208 - 


3.6625 -6.7959 1.5681 2.5541 


14.716 27.443 


60 -54.117 -80.899 


-9.9713 -13.594 


-35.514 -53.091 


-4.2542 -5.8025 - 


4.5432 -6.7939 1.9029 2.5918 


18.356 27.441 


90 -67.458 -80.910 


-9.1111 -13.607 


-44.267 -53.098 


-3.8893 -5.8080 - 


5.6594 -6.7948 1.7357 2.5944 


22.885 27.445 



Z 




FIG. 10. Plot of the mapping function z(y) defined in Eq. | |B5[ |. 
The four curves, each with 90 Gauss points, correspond to bo = 
10 (small solid triangles), bp = 1-5 (open boxes), bg = 1-2 (solid 
boxes), and bo = 1.1 (small solid circles). 



NN fits, the meson propagators have the form 



Vb{z) = F{z) 



mt 



il + A\zQ~z\ 



K 



n4 



Kl + A\zo-z\ 



(B3) 



where mi, and A^ are the meson mass and form factor 
mass of meson b, respectively, and 



A = 2pp' 



zq 



^:!^i^^...,, ,B, 



with ^ = if the incoming and outgoing momenta (p 
and p') are equal in magnitude, and ^ > when p ^ 
p' . [Since zq > 1, the absolute value in (B3) can be 
ignored.] At the largest momenta used in the numerical 
solutions, A ~ 1800 GeV'^, compared to (for example) the 
pion mass squared m^ ~ 2 x 10^'* GeV^. Under these 
conditions, numerical evaluation of the integral, with its 
narrow and large forward peak is difficult. 

To carry out the integration numerically, we apply a 
Gauss-Legendre quadrature rule to discretize the inte- 
gration variable z into a grid of Gauss points distributed 
between -1 and +1. However, it is convenient to map 
the integration variable z ^- y, where y will vary slowly 
in the region where F{z) is large. One possible func- 
tion that will accomplish this, and still remain bounded 
between —1 and 1 is 



z{y) 



ta.nh[y/{bo- 1)] 
tanh [l/(6o - 1)] 



(B5) 



When bo is close to 1, this function spreads out the in- 
tervals near z = ±1, so that more Gauss points may 
be placed in the region of the forward peak. The map- 
ping function, for four different values of bo, is shown in 
Fig. 10 Using these mappings, the propagator becomes 



a function of the integration variable y (instead of the 
physical variable z), and the integrand, as a function of 
the mapped variable y, is shown in Fig. ni] (for the same 
four choices of b) . The angular integral is transformed to 



' r' dz 

dzF{z)= / dy-F[z{y)]. (B6) 

1 j-i "y 



Clearly, without any mapping (6o > 10) the forward 
peak is very poorly sampled. If we take b^ = 1 + e, with e 
small, but larger than zero (to avoid the "singularity" at 
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FIG. 11. Plot of the meson propagator F (with m^ = 0.14, A;, = 4, zq = 1, and A = 1800) as a function of the mapped variable y defined 
in Eq. | |B5| . The four curves, each with 90 Gauss points, correspond to bo = 10 (small solid triangles), feo = 1-5 (open boxes), feo = 1-2 
(solid boxes), and bg = 1.1 (small solid circles). The right-hand panel covers the region shown in the upper right corner of the left-hand 
panel. 




FIG. 12. Plot of the mapping function z{y) defined in Eq. ( |B8| , 
with A = A. The three curves, each with 90 Gauss points, corre- 
spond to zo = —0.5, zq = 0, and zq = 0.5. Each curve fiattens out 
in the region of z = j; = zq. 



&0 = 1), we are able to evaluate the integrals accurately. 
In practice, we found that 6o = 1-2 (one of the cases 
shown in the figures) gave a good sampling of the forward 
peak and the integration converged rapidly. 



3. Mapping for the exchange terms 

The exchange propagator has the same form as 
Eq. (B3), except now 



zq 



p^ + p'^ - {Ep + Ep, - wy 



2pp' 



-^ -1 



2W 
P 




FIG. 13. (Color online) Plot of the meson propagator F (with 
mi, = 0.14, At, = 4, and A = 1800) as a function of both the 
unmapped variable y and the mapped variable z{y), defined in Eq. 
| |B8| , for each of the choices zq = —0.5, zq = 0, and zq = 0.5. The 
90 Gauss points distributed between —l<y<l are shown on each 
curve. All the curves peak at y = zq, with the unmapped functions 
F{y) having narrower peaks than the mapped function F{z(y)). 



where W is the total energy of the NN system in its 
center of mass, and the last expression holds if p = p' — > 
c». We see that it is now possible for j^ol to be less than 
unity, so that the peak can lie in the region of integration. 
To handle this case we need a mapping where the num- 
ber of Gauss points is strongly distributed near z = zq, 
where F{z) is sharply peaked. It is convenient to use 
a mapping which maps the points y ~ zo,l,— 1 into 
z — zq,1,—1 respectively. Such a mapping is 



(B7) 



ziy) 



z+iy,b) Zo<y<l 
z-{y,b) ~1 <y < zq 



(B8) 
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where 



^ tanh[C(2^-l-zo))/(l-zo)] 
z^iy,zo)^,{l~zo) ^^^^ 

+ 5(^0 + 1) 
z (v z-)-U\ \ z\ ^^^^^^(2^ + 1 - zo))/(l + zo)] 



+5(^0-1) 



(B9) 



and C is a parameter. This mapping (with (^ = 4) is 
shown in Fig. 12 Note that each curve distributes the 



90 Gauss points so that they are densely spaced in the 
region of z = Zq. 

Both the mapped and unmapped meson propagators 
are shown in Fig. [iSJ for each of the cases zq = —0.5, 
0, and 0.5. Note that the unmapped propagator has a 
very narrow peak at y = zg, and that this peak is de- 
scribed by only a few Gauss points, even when the total 
number of points is quite large (90 in this case). The 
mapped functions, however, have their peaks spread out 
with many Gauss points distributed around the peak. 
From this picture, it is easy to see why the mapping 
allows a more accurate evaluation of the integrals. (In 
practice, a better result was obtained by dividing the 
interval into (.1 = [— l,Zo] and (.2 = [zo,l], and doing 
separate Gaussian integrations over each interval. The 
points assigned to each interval, n^, depended on the 
value of zq. If zq < —0.5, we chose rii — n/4 and n2 = n, 
for —0.5 < Zq < 0, ni = n/2 and n2 = n — n/A, for 
< Zq < 0.5, ni = n — n/A and 712 = rt/2, and for 
0.5 < Zq, Til = n and 71,2 = n/A, so that in all cases we 



actually used 5n/4 Gauss points.) 

Finally, the angular convergence is illustrated in Tables 



XV and 



kernel v'l 



XVI for a typical matr ix element [the isoscalar 
^ defined in Eq. (A12|] at one of the most dif- 
ficult points: the diagonal elements with p = p' — Pmax 
where Pmax — 30 GeV, the largest momenta used in the 
numerical solutions of the equations. Table |XV| shows 
the convergence of the direct terms (those with xq = 1) 
and Table XVI the exchange terms (those with xo = —1). 
Note that the direct terms calculated with unmapped in- 
tegrations (the columns labeled with y in the tables) do 
not converge, and that the convergence of the exchange 
terms is very marginal. With the mappings (columns la- 
beled z), the direct terms converge even at n = 20, but 
excellent results for the exchange terms requires n > 40. 
Because the high-momentum matrix elements only make 
a small contribution to the solutions of the equation, 
smaller values of n give very good results for the phase 
shifts, and these improved angular integrations only af- 
fect the accuracy with which the asymptotic deuteron 
wave functions can be determined. 



Appendix C: DEUTERON WAVE FUNCTIONS: 
SOME DETAILS 



In this appendix we start by showing in detail how 
the matrix equation (3.26) is derived, and then use the 



equation to study the asymptotic behavior of the wave 
functions, obtaining the results (3.41 ). We conclude with 



a demonstration of how the general normalization condi- 
tion (3.12) reduces to (3.39). 



1. Equations for the partial wave deuteron w^ave functions 



The equation for the deuteron wave functions, in partial wave form, can be extracted from the general result (E36) 
of Ref. I. This equation uses a kernel antisymmetrized under particle exchange (with the superscript J ~ \ suppressed) 



vxtt'kx'Sp^ k) - \[v:,:i'X:i,,,{^.Po, k) + s^v2t:t:k^^(p^ -p^-^ k)}. 



(Cl) 



We will show how the full kernel with well defined parity, Eq. ( Al I, will emerge automatically, but the phase of the 



amplitude under particle interchange cannot emerge automatically in the CST (as discussed in Refs. I and II); that 



phase must be imposed by hand from the start by using the symmetrized kernel (Cl I. The phase Sp^ = ±1 is related 



but not equal to the phase under particle exchange [because the hehcities Ai, A2 are not exchanged in (Cl)] 
Using this kernel, the deuteron bound state equation is 



m 



^\,L,xAp) 



k^dk MT^ —+p2^+p 



(27r)3 



E KZ 



x.,,.,Mk)G''ik)-T;^^^^,^{k) 



(C2) 



The dependence of the amplitudes on the total momentum P has been dropped, and the energy factors of m/E, 
written explicitly, will be absorbed later when we go to the matrix notation of Eq. (3.26). At this stage the sum is 
over all helicities /ii , ^2 and the intermediate rho-spin p. 

The symmetries of the vertex function are a direct consequence of the symmetries of the kernel. Under parity, the 
deuteron kernel transforms as 



rv. 



+ P2: + P 

d, AiA2,piM2 



(p, k)r- 



PP2V d!'W^,-p^^p^,_{p, k). 



(C3) 
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Hence the equation satisfied by the transformed vertex function, 



"^ r+P2 / \ I "~ k dk 



(2^) 

k'^dk 

(2^ 



Y.K'X->.......Mmk)^Ti'^ 



M1M2P 



m 



Ek 
+p 



(C4) 



is identical to the original equation, showing that the parity relation (3.24) is satisfied. Using this result we may 



reduce the sum over 1x1,^2 from four to two terms. Dropping the redundant pi — p'x = + indices, we can write the 
summation over the two values of pi explicitly and use the parity relations to reduce the equation to a form in which 
only /ii = +\ appears: 



-1P2 



°° k^dk m 
j2^E~k 






P2P 



l^HLyYP-f (<5„„,p)G"'(fc)r% , (fc) 

^ ^ P2P 



P2:>^d 



(fc) 



(C5) 



where in the last line we have recovered the fully symmetrized kernel of Eq. (Al). The deuteron channel must have 



the phases SpQ = P2, appropriate to an isospin state, and Ss = P2 = P- Using the notation V = V^^j^ g (and, for 



the helicities, ±1 — ^ ±) leads to the identifications 



Vl\+4p,k)=vl^ + 



-+ 



V 
V 



P2 - 



P2 - 



_{p,k) = v^ 
(p, k) = v^ 



+ - 



Y'^^X^P, k) = vll + Y'[^-+{p, k) = vl^ - 



VZ+_(p,fc)=«r6+ VZ+^_(p,fc)=z;4-- VI+_(p,fc)=i;r4+ VZ-_(p, fc) - ^3- - (C6) 

, |r7_), respectively. Note that they 



where the kernels are written in the order for coupling to |r^_|_), |r^_|_), |r 



reproduce the matrix (3.28), finishing this discussion. 



2. Asymptotic behavior of the wave functions 



Starting from Eq. (C5), the wave function (3.22) has the following general behavior at large p 



r. 



m . 

^P 



x.X2-xM-^'iG''{p)j^Kx,^xAp)^G''{p) I dkV,{p,k)C{k), 



(C7) 



where, in this section, we return to the notation k = |k| 
for the magnitude of the three-momentum, and V^ is 
the appropriate partial wave projection of the Feynman 
OBE amplitudes [shown in detail in Eq. (C6)]. All we 
need know about the function C(fc) is that it provides 
the convergence for the k integral. The kernels include a 
factor of m/Ep, which contributes to the large p behavior. 



The kernels Vd are multiplied by nucleon form factors 
[recall Eq. ( |3.5| )], one of which depends on the intermedi- 
ate momentum (fc) and the other on the final momentum 
(p). The final state form factor, h{p), plays an important 
role in the asymptotic behavior of the wave function. It 
depends only on p'^, which is the square of the mass of 



the off-shell nucleon 



r = {Md 
p- ^ 



Epf 



^ ^ 2Mdp for p > m . 



Hence, 



Kp") 



(A 



,2\4 



16M4 p 



4,^4 



for p 3> m. 



(C8) 



(C9) 



so the form factor behaves like p~^. This contributes a 
k^^ behavior to the function C(fc), which, together with 
the internal nucleon form factor h{k) that is part of V^;, 
assures that the integral over k will converge. Hence we 
may conclude that the wave functions go like 



V'AiA2,A,(p) 



GP{p)Ma{p.~k)C{k). 



(CIO) 
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where the mean value theorem has been used to write the 
integral over k, with k a four-vector whose three-vector 
part has some fixed (unknown) value k in the interval 
[0,oo). 

To finish this discussion we examine the asymptotic 
behavior of the partial wave amplitudes for one of the 
OBE terms in the kernel. If the form factors h{p) and 
h{k) are removed, this integral has the generic form 



Dbip,^) = 



\ , . m f{Ab,q) , . 



(Cll) 



where q^ is the momentum transferred by the boson with 
mass TTif), form factor /(Ab, g), and momentum dependent 
coupling (7b(p)— > p^ (as p— >■ oo, where I is yet to be deter- 
mined). We have isolated the m/Ep term incorporated 
into the definition of Vj, and d(z) includes the z depen- 
dence that arises from the generic rotation matrices of 



the type tabulated in Table XI In Ref. I the meson form 
factors have the form 



/(Ab,<z) = 



A? 



A? 



(C12) 



where rif, = 4 for all mesons. 



The asymptotic behavior of the integral ( Cll ) depends 



on whether or not it is a direct or an exchange term. The 
momentum transferred by the mesons in the direct terms 

q^^ {Ep ~ Ekf - (p - k)2 

-^ 2m^ - 2p{Ek - kz) (as p ->cx), k fixed) (C13) 

and is never positive, so the integrand is a smooth func- 
tion of z. The asymptotic behavior, for a fixed k^kg, is 
therefore 



A'"(p,k) 



,"6-1-2 



p 



(C14) 



As shown below, these terms are vanishingly small 
compared to the exchange terms. 

The exchange terms have a very different structure. 
Here the momentum transferred by the meson is 

q^ = (Md-Ep-Ekf-ip + kf 

^ m2 + (Ma - Ekf - k^ - 2p (M^ - Ek + kz) 



— 2p(Afrf — Ef; + kz) as p — > oo . 



(C15) 



When p is very large, this is zero along a line in the k, z 
plane given approximately by Md — Ek+^z = 0. Along 
this line the integrand has a cusp which sits at 



Md - Ek 



(C16) 



in the region of integration (— 1 < z < 1) for values of k 
in the interval 



Ml 



2Md 



< k < oo . 



(C17) 



Contributions from the cusp clearly dominate the inte- 
gral, because in that region the integrand is constant. 



Hence, letting rrib — At, in order to simplify the deter- 
mination of the asymptotic p-dependence, the exchange 
integral is approximately 



i?-(p,k) 



dz 



dz 



A? 



Ag 



"6 + 1 



d{z) 



A? 



A2 + 2pk(z, - z) 



"b+i 



dz 



9A2 



A2 
^^6 



A? 



2pk(z — Zc 



d{z) 

"6 + 1 



ribpk 



d{z) 
(C18) 



where we have integrated by parts, keeping only the lead- 
ing term (the corrections go like higher powers of 1/p). 

Collecting these results together we see that the re- 
duced partial wave amplitude, Db, behaves like 



^(p,ko 



P 



t-{2+nt) 



direct 
exchange 



(C19) 



where nb is the power of the form factor at the NNb 
vertex and i the power of the p dependence of the NNb 
vertex function. 

It remains to determine the power L First, note that 
the s ingle particle p-spin matrix elements given in Table 



XIV 



each contribute a factor of \/j>/{2ra) [from the high 
momentum limits (A5|]; their product gives one power 
of p implicitly contained in Eqs. (A45)-(A48|. 



Next, examine any additional p dependence coming 
from the two particle p-spin matrix elements given in 
Eqs. (A45|-(A48). There are terms of order p coming 



from the off-shell energy factors A_Ej defined in (A24|, 



but the condition (A25l shows that these factors cannot 
be multiplied together to make a term of order p^. Spe- 
cial terms of order p^ seem to be present in the vector 
meson exchange terms, but we will show below that these 
terms cancel. Hence the final result is the product of two 
terms of order p giving 1 = 2. 

To show the cancellation of the special p^ terms con- 
tributing to vector-meson exchange, consider the three 
terms in Eq. (A47) which give rise to them: the prod- 



uct of Rvi's (first term), the p^ multipying the D4 terms 
(fourth term), and the last term. It is sufficient to look 
at the exchange-term matrix elements where particle 1 is 
off-shell in the final state, so that AE2 = 0. At large p 
use the fact that AEi — )■ (pi + l)p to obtain 



RviRvi — ^ P 



2 9v'^v 



(2to)2 

2 ( 9v>^v\ 



2 



111 



l-'-2 



I1I2 



RvsDa -^ p 



2 9v^v 

(2m)2 



{n,-v„){pi + l)2\^Tll2. (C20) 



Adding these together, and replacing pi by p2 [in £^c- 
cordance with the notation for the exchange terms in 
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FIG. 14. The p dependence (for k fixed at 816 MeV) of the 16 model WJC-1 deuteron kernels Vi(p, k) defined in Eq. | |C6| (but with the 
nucleon form factor s set to unity). The left panel shows clearly that the total kernel (direct plus exchange) approaches a constant at large 
p, as predicted by | |C19| with £ = 2. The left panel shows that the asymptotic p dependence of 5 of the direct kernels behaves like p"**, 
again predicted by ||C19| with nj, = 4 and £ = 2. The other direct kernels decrease more rapidly, presumably because they have an £ < 2. 



Eq. (A12|] gives the following compact result 



3. Reduction of Eq. (3.12) to Eq. (3.39) 



gv{p) ^ p^ 77t::^i'^v - v^){p2 + 1) 



(2m)2 



li + 2\itI 1 



(C21) 



Note that this is zero if p2 = ^1, insuring immediately 



that 8 of the deuteron kernels ( C6 1 have no special terms 
of order p^. To see that this is also true of the other 8 
cases set /O2 = 1 (and recall that p'^ = 1) and substitute 
for th e ma trix elements 1 and r^ using the results from 
Table XIV (with p2 ^ Pi) and the high momentum limits 
dASl) 



yviv) ^ 2p^ 



QyKy 



2m {2111)^^"^ ''''^ 
x(l-2Ai)(c'-.s')l2. 



(C22) 



^, which is true of all of the am- 

We have 



This is zero if Ai 

plitudes (and, in particular, the remaining 
proved that the special p^ terms coming from vector me- 
son exchange cancel, and that 1 = 2. The conclusion is 
that the wave functions should go like 



zP{v)^G''{p)h{p')D^-{^,k^) 



(C23) 



which gives the result (3.411 as expected 



To confirm the predictions of Eq. \C&\ , Fig. Il4j shows 
how the 16 deuteron kernels V^i(p, k) vary with p for a 
fixed k . 



Substituting the expansion (3.16) into Eq. (3.12), and 



using the orthogonality of the Dirac spinors, compactly 
written as 

uP{-^, X)fu''' (-P, A')= ^SyxSp'p (C24) 



gives 



/,pt 



5^,y^= /d^p<;,^,,(p)^^^,^^,Jp) 






where the V 



X^^iA,,A,(P'), 
T+P.+P' 



(C25) 



V 



includes the factors of 



m^/EpEp, [cf. Eq. (E13) of Ref. I]. 

Next we express the helicity amplitudes in terms of 
the wave functions Z£{p) using (3.30). To this end recall 



Eq. ( A6 ) for the two-component helicity spinors and use 
the convenient relation icr2Xx — ~2Ax_^ and the matrix 
elements 



-xl 



X^^^cr • pio-2X;,^= Xl;,^X-Ai = '^Ai,A2 



(C26) 
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where A = Ai — A2 and the rotation matrices are 

4q(6() =cos0 

dloi&) = dl-M = -dlM - -dUfiid) = -75 sin0 

d\M = dl^^_^{e) = \{i + cose) 

dl,^M^dl_M = W-^ose). (C27) 

The result is 



tion condition (3.201 and the definitions (3.211 



V'++,A.(P) = y/ij_^dz\dl^„m^ {u + V2w) 



\/6 



1 '- ' V2v 



U + V zw 



V^+-,A.(P) = V^/^ dz ML,i Wl' (^^" - ^) 



\/6 



\/2u-^ 



^^^^^-^^^P^ = ^ H-TI'^^O ^'"'^-^^'^ V';+,..(P) = -^ frf^MLoWP-.^-^ 



+^....^<.W 



^+-aM--7^s dz\dl^,{e)\'v, = -j^^ 



Vt . 



V'A.A.A.lP) = -2An/- <5a,.a. «.(p) ~ ^ rfL.o(^) 



t^t(p)A ,1 



V2 



(C31) 



+^^/2'^'4.aW 



(C28) 



ExpHcitly, for each individual case, using A = ±i 
^Ia,a,(p)= ^ ("(P) + V2w(p)) dl^,{e) 

V'aVa,a.(p)= ^ ("(p) - ;7I^(P)) '^I-^aI^) 



4. Coordinate-space expansion functions 

The coordinate-space wave functions discussed in 
Sec. |III G 3| are constructed from the spherical Bessel 
transforms (3.37), which can be written in terms of the 



Jo Bessel transform 



fe.A.(p)=-2AJ^i^.(p)dL,oW 



G\{r) = rJi p'dpu{pr)G\{p) 
Jo 

'^^l+l ( __\ I „2 

r dr 



V'A,-A,A,(P)= -2AJ — l.t(p)dL,2AW- (C29) 



We now can reduce the first term in the normalization 



y+^ (\ d 



P dpjQ(pr) ^- 

P 



(-1) 



Af, \r dr 



Ghir) 



(C32) 



condition ( C25 1 . Takin g A^ = A^ and using the normal- 
ization condition (3.20) for the d^ functions gives 



Hence, for S and P states it is sufficient to calculate the 
following Fourier sine transforms: 



Ghir) 



pdpsin{pr)Gl,{p) 



AAe-'^-e-^'[l + ^Z,{l-R^)] 



1 =27r / p'dp / sinM42|V,["A.A,(p)l' + 2|^|._,,A.(P)I 
+2|^aa,a,(p)I'+2|^Ia,a,(p)I'} 
= -y p^dp\^[u{^)^^/2w{v)) 

+2 {u{x,) - -^ w(p)^ + Zvl{v) + 3^;2(p)} 

p^dp\%?{^) + ii;2(p) + w2(p) + i;2(p)| (C30) where z,, Z,, i?,, and A, were given in Eq. ([3^ and the 



Goir) = ^Ji pdpsm{pr)G^{p) 
Jo 

= ^M^Z^K,{Z,,) 



(C33) 



in agreement with Eq. (3.391. The derivative terms can 



be similarly reduced and expressed in terms of the kernels 
given in Eq. ( C6 ) . 



modified Bessel functions of the second kind are 



Kn{z) = j^ ^ dte-''\t^~iy'--. (C34) 



Before leaving this section we note that the partial- 
wave deuteron wave functions follow from the normaliza- 



(2n-l)!!yi 



The P-state wave functions are obtained by diffcrenti- 



ation, as outlined in Eq. (C32). They are 
-Gl{r) 



G\{r) 



r d 
W,dr 



AAR,e- 



1 



-,-Zi 



1 



2^4 (l - i?j) 



G«(r) 



r d \G^{r) 



Stt 



M^Z'„KoiZ,,). (C35) 



first transforming the functions 

gUp) -- 
G'SAp) - 



Finally, the D-state wave functions are computed by which gives 



m; 



MP 



^(m|+p2)(M2+p2)3 
2 Ml 



7r(M2+p2)7/2 
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(C36) 



G'oAr) 
GUr) 



f J pdpsm{pr)G},M = B.l e^^' - e'^^ - |Z,e-^' 



2\^ 



(5 - 6i?2 + Rj) +Z,{1~ i?2) 



2 / pdpsin(pr)GS<i(p) = ^^Af2zf,i^2(^n 





IStt 



where Bi was given in Eq. (3.52). Differentiating these following (C32) gives 

2 ( 



^^W-i^(d^!:) Go.M = ^^^^?e- 



1 



^r 





-e-^ 




2 


Id 11 


2 


^^ 


(ir r 





i + |- + ^ + i(i-i?D(i + ^0 + i^f(i-i??)' 



G?M 



At large r the functions G^ go like 



G'iAr) 



Ibn 



M,iZ'^Ko{Zn). 



r 



plied by a fractional power of r 



liniG[,Hr) = ^y^Af^Z3/2e-^" 

r-foo 6 \ t: 

\\n.GUr)^\J^-MlZl'-e'^- 

r-i-oo 3 V TT 

limG^W^iy^M^Zf/^e-^" 

r->oo 15 V TT 



(C37) 



(C38) 



(C40) 



lim Gl{r) = A, e"^' 
lim G\{r) = A,Rie-''' 

r—^oo 

lim Gl{r)^B,Rfe''^ 



At small r the functions G^ have the expected r^+^ 
Zf+i behavior 







[i + l 


3l 

+ 72 



lim G^(r) = iv4,Z, 1 - i?. 



i-->-0 



(C39) 



lim Gl(r) = iA,Z2(l - i?,)2(l + 2i?,) 

lim G^(r) = 4L5,^^3(^ _ ^^)3(i ^ 3^^ ^ 8^2) (^41) 



r-s-O 



However, the tail functions contribute some nonanalytic 
behavior at small r: 

limG^(r) = ^M2z„. 



T-->-0 



37r 



This is the asymptotic behavior expected for solutions of 
the Schrodinger equation with orbital angular momen- 
tum £. The tail functions fall-off like exponentials multi- 



lim Gr ir)^- — M^Z'J log Z„ - log 2 + 7 

lim G^ (r) = - ^M^Z^ ( log Z„ - log 2 + 7 ) , (C42) 
where 7 — 0.5772 is Euler's constant. 
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